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Abstract 



We make a high-precision Monte Carlo study of two- and three-dimensional 
self-avoiding walks (SAWs) of length up to 80000 steps, using the pivot algo- 
rithm and the Karp-Luby algorithm. We study the critical exponents v and 
2A4 — 7 as well as several universal amplitude ratios; in particular, we make 
an extremely sensitive test of the hyperscaling relation dv = 2A4 — 7. In two 
dimensions, we confirm the predicted exponent z/ = 3/4 and the hyperscal- 
ing relation; we estimate the universal ratios {R"^) / {Rl) = 0.14026 ± 0.00007, 
{R'i)/{Rl) = 0.43961 ± 0.00034 and = 0.66296 ± 0.00043 (68% confidence 
limits). In three dimensions, we estimate z/ = 0.5877±0.0006 with a correction- 
to-scaling exponent Ai = 0.56 ± 0.03 (subjective 68% confidence limits). This 
value for z/ agrees excellently with the field-theoretic renormalization-group 
prediction, but there is some discrepancy for Ai. Earlier Monte Carlo esti- 
mates of z/, which were ~ 0.592, are now seen to be biased by corrections to 
scaling. We estimate the universal ratios = 0.1599 ± 0.0002 and 

= 0.2471 ±0.0003; since > 0, hyperscaling holds. The approach to is 
from above, contrary to the prediction of the two-parameter renormalization- 
group theory. We critically reexamine this theory, and explain where the error 
lies. 
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1 Introduction 

The self-avoiding walk (SAW) is a well-known lattice model of a polymer molecule 
in a good solvent [1]. Its equivalence to the n = limit of the n-vector model 
[1,2,3,4,5,6,7] has also made it an important test-case in the theory of critical phe- 
nomena. 

In this paper we report the results of an extensive Monte Carlo study of two- 
and three-dimensional SAWs of length up to 80000 steps, using the pivot algorithm 
[8,9,10].^ We make a high-precision determination of the critical exponents v and 
2A4 — 7 as well as several universal amplitude ratios. In particular, we make an 
extremely sensitive test of the hyperscaling relation d.v = 2A4 — 7, which plays a 
central role in the general theory of critical phenomena (Section 1.1). 

Our results have also led us to reexamine critically the conventional theory 
of polymer molecules, the so-called "two-parameter renormalization-group theory" 
[11,12]. Indeed, such a reexamination is unavoidable, as our Monte Carlo data are 
inconsistent with this theory as it has been heretofore applied. But this is because 
— as we explain in Section 1.2 — the theory has heretofore been applied incorrectly! 



^These computations were carried out over a 4-year period on a variety of RISC workstations. 
The total CPU time was several years, but we have by now lost track of exactly how many! 
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These points were first made three years ago by Nickel [13], in an important but 
apparently under-appreciated paper; they have recently been extended by one of us 
[14,15]. 

1.1 The Problem of Hyperscaling 

One of the key unsolved problems in the theory of critical phenomena has been 
the status of the so-called hyperscaling relations (scaling laws in which the spatial 
dimension d appears explicitly). These relations have long been known to rest on a 
much more tenuous physical basis than the other scaling laws [16,17,18,19,20,21,22]. 
Indeed, it has been understood since the early 1970's that hyperscaling should not 
hold for systems above their upper critical dimension c/„: for d > du the critical ex- 
ponents are expected to be those of mean-field theory, and these exponents satisfy 
the hyperscaling relations only at d = du-^ (For models in an n-vector univer- 
sality class, including the SAW, c/„ equals 4.) It has generally been believed that 
hyperscaling should hold in dimensions d < du, but in our opinion there is no partic- 
ularly compelling justification for such a belief (although the claim itself is probably 
correct). Hyperscaling in the three-dimensional Ising model is the subject of a con- 
troversy that has been raging for 25 years, and which is still not completely settled."^ 
We remark that hyperscaling is also of interest in quantum field theory, where it 
is equivalent to the nontriviality of the continuum limit for a strongly-coupled (^'^ 
field theory [48,49,50,37,51,52,7]. 

Although the hyperscaling relations appear naively to be ineluctable conse- 
quences of the renormalization-group approach to critical phenomena, closer exami- 
nation reveals a mechanism by which hyperscaling can fail: the so-called "dangerous 
irrelevant variables" But the much more difficult question of whether this viola- 
tion actually occurs in a given model can be resolved only by detailed calculation. 
Unfortunately, a direct analytical test of hyperscaling appears to be possible only 
at or in the immediate neighborhood of a Gaussian fixed point, that is, for asymp- 
totically free theories [59,60,61,62] or for small e = du — d or large n [63,64,65,66]. 
We note that the real-space RG [67,68] and field-theoretic RG [69] frameworks, 
as typically used in approximate calculations, implicitly assume the hyperscaling 
relations, so they cannot be used to test hyperscaling. 

It is therefore of interest to make an unbiased numerical test of hyperscaling, 
working directly from first principles. One approach is series extrapolation [70,71], 
which affords a direct test of universality and scaling laws, including hyperscaling. 

^This belief has now been confirmed by rigorous proofs of the failure of hyperscaling for the Ising 
model in dimension c? > 4 [7], the self-avoiding walk in dimension c? > 5 [1,23,24], and spread-out 
percolation in dimension c? > 6 [25,26]. 

^See [17,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42] for series-extrapolation work, and 
[43,44,45,46,47] for Monte Carlo work. 

■*This mechanism was proposed independently by Fisher [53] and Wegner and Riedel [54, p. 250 
and footnote 8] in the early 1970's. For further discussion, see also Ma [55, Section Vn.4], Amit and 
Peliti [56], Fisher [57, Appendix D] and van Enter, Fernandez and Sokal [58, Section 5.2]. 
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It gives numerical results of apparently very high accuracy: the claimed (subjec- 
tive) error bars on critical exponents are on the order of ±0.001-0.005 [71], which 
is comparable to the best alternative calculational schemes. However, as is inherent 
in any extrapolation method, the results obtained depend critically on the assump- 
tions made about the singularity structure of the exact function, notably the nature 
of the confluent singularity, if any [33,38,41,72,73,71]. Indeed, estimates by differ- 
ent methods from the same series sometimes differ among themselves by several 
times their claimed error bars. This, together with systematic differences between 
lattices of the same dimension, accounts for much of the controversy over hyper- 
scaling. Quite a few extra terms would be needed to resolve these discrepancies 
in a convincing manner [33]. Unfortunately, the computer time required to eval- 
uate the series coefficients grows exponentially with the number of terms desired, 
while the extrapolation error is proportional to some inverse power of the number 
of terms (the power depends on the details of the correction-to-scaling terms and 
the extrapolation method). 

In a Monte Carlo study, by contrast, one aims to probe directly the regime 
where the correlation length is ^ 1. (For SAWs this corresponds to a chain 
length iV ^ 1.) The method affords a direct test of universality and scaling laws, 
including hyperscaling. In practice, however, it has been extremely difficult to 
obtain good data in the neighborhood of the critical point. There are two essential 
difficulties: finite system size [74,75,76] and critical slowing-down [77,78,79]. For 
spin models and lattice field theories, these two factors together imply that the CPU 
time needed to obtain one "effectively independent" sample grows as ^L'^^^ > i^''"''^, 
where d is the spatial dimension of the system and z is the dynamic critical exponent 
of the Monte Carlo algorithm.^ (This situation may be alleviated somewhat by 
a new finite-size-scaling technique [80] that yields accurate estimates of infinite- 
volume quantities from Monte Carlo data on lattice sizes L <^ i^.) The situation 
for SAWs is rather more favorable: one can simulate a SAW directly in infinite 
space, with no finite-size corrections or L'^ factor in the CPU time. There is, to be 
sure, critical slowing-down; but vast progress has been made over the last decade 
or so in inventing new and more efficient algorithms for simulating the SAW [81]. 
In particular, using the pivot algorithm [8,9,10] one can generate an "effectively 
independent" iV-step SAW (at least as regards global observables) in a CPU time 
of order ~ iV [10,81]. (This is the best possible order of magnitude, since it takes 
a time of order N merely to write down an iV-step walk!) Since ~ iV, this 
corresponds to a CPU time ~ ^^^'^ in a spin system — which (if d > 2) is better 
than ~ 1^''"'"^ even if 2; = 0! So the SAW is a uniquely favorable "laboratory" for 
studying the problem of hyperscaling. 

^Conventional local algorithms have z 2, while the new collective-mode algorithms [77,78,79] 
can have z <^ 2 and in some cases even z = 0. 
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1.2 Which Quantities are Universal? 

Over the past four decades, various mathematical models have been employed 
to describe the behavior of linear polymer molecules in a good solvent.^ Among 
these models are the self-avoiding walk [1], the bead-rod model [82], and the con- 
tinuum Edwards model [83,84,85,86,87,11,12]. The detailed behavior depends on 
the specific model chosen, just as the detailed behavior of real polymer molecules 
depends on the particular chemical structure of the polymer and solvent (and on the 
temperature). However, it has long been understood that some aspects of polymer 
behavior become universal in the long-chain limit iV — )• oo (where N is the number 
of monomers in the chain). Unfortunately, there has been considerable confusion 
about which quantities are universal and which are not. In this subsection we sum- 
marize recent work of Nickel [13] and one of us [14,15] which clarifies this issue. 
(For further discussion, see Section 5.2 below.) 

Standard renormalization-group (RG) arguments predict [13] that the mean- 
square end-to-end distance (-Rg), the mean-square radius of gyration (-R^) and the 

second virial coefficient A^""'^ = {N^ M^^^^^^J NAvogadro)A2 of any real or model 
polymer chain should have the asymptotic behavior 

(Rl) = An^N'^{l + b^^lN-^^+...) (1.1) 

{R') = A^^iV^ll + 4>-^^ + • • •) (1-2) 
4"°') = AAN^'^il + b^j^N-^' + ...) (1.3) 

as iV — > oo, where d is the spatial dimension.'' The critical exponents v and Ai 
are universal. The amplitudes A^^^, A^^^, A^, 6^], are nonuniversal; in fact, 

even the signs of the correction-to-scaling amplitudes 6^], [and their various 

combinations such as b^^^ = — (c//2)6^^^] are nonuniversal. However, the RG 

theory also predicts that the dimensionless amplitude ratios Ar^/Ar^^ Aa/A'K^, 

^rI/^rI /^rI universal [13,88]. 

So there is no reason why the correction-to-scaling amplitudes should have any 
particular sign. In the continuum Edwards model, the effective exponents Veff,Re = 
I (/ log (i?g)/c/ log iV and Veff,Rg = \d\og{R^g) / d\og N and the interpenetration ratio 
\E' = 2(d/127rY^^A2""^'^ / {R'^)'^^^ all approach their asymptotic values from below 

[11,12,89,90,91]: that is, > and < 0. On the other hand, high- 

precision Monte Carlo data on lattice self-avoiding walks — see Section 4 below, 
as well as Ref. [13] — show clearly that these quantities approach their asymptotic 
values from above; and the same occurs in the bead-rod model with sufficiently large 
bead diameter [92]. Indeed, this behavior is almost obvious qualitatively: short self- 
avoiding walks behave roughly like hard spheres; only at larger N does one see the 

®Here "good solvent" means that we work at any fixed temperature strictly above the theta 
temperature for the given polymer-solvent pair. 

^In (1.3) we have assumed for simplicity that the hyperscaling relation di^ = 2A4 — 7 is valid. 
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softer excluded volume (smaller characteristic of a fractal object. In any case, 
all these models are in excellent agreement for the leading universal quantities i^, 
ArJAr^ and ^* = 2{d/127rY/^AA/A'j(\ and they are in rough agreement for the 

universal correction-to-scaling quantities Ai, h^R^/h^^^ and I^^Re- 

It is thus misguided to analyze the experimental data in the good-solvent regime 
by attempting to match the real polymer molecules to the continuum Edwards 
model via the correspondence ZEdwards = o,N^^'^ (where a is an empirically deter- 
mined scale factor depending on the polymer, solvent and temperature)®: the con- 
tinuum Edwards model can predict only the universal quantities. Indeed, there is 
evidence [93,94,95,96] that real polymers in a sufficiently good solvent behave like 
self-avoiding walks, i.e. they approach \E'* from above; in this case they cannot be 
matched to any value of ZEdwards- This behavior has heretofore been considered 
paradoxical; in fact, it is quite natural. (Huber and Stockmayer [94] attributed this 
behavior to the effects of chain stiffness. In fact, as pointed out by Nickel [13], chain 
stiffness is quite irrelevant here, as the effect occurs also for perfectly flexible chains, 
such as self- avoiding walks or the bead-rod model.) 

These points have been made previously by Nickel [13]. Similar comments have 
been made with regard to liquid-gas critical points by Liu and Fisher [97]. 

In summary, the error of all two-parameter theories is to fail to distinguish 
correctly which quantities are universal and which are non-universal. In particular, 
the modern two-parameter theory begins from one special model — the continuum 
Edwards model — and assumes (incorrectly) that it can describe certain aspects of 
polymer behavior (e.g. the sign of approach to \E'*) which in reality are non-universal. 

Remark. A very different limiting behavior is obtained (in dimension d < 4) 
if we take simultaneously iV — )• oo and T —i- Tg such that x = N'^{T — Tg) remains 
fixed, for a suitable crossover exponent (p. In a separate work [14,15], one of us has 
argued that it is precisely this universal crossover scaling behavior in an infinitesimal 
region just above the theta temperature that is described by the continuum Edwards 
model. 

1.3 Plan of this Paper 

The plan of this paper is as follows: In Section 2 we review the needed back- 
ground information about the self-avoiding walk and the pivot algorithm. In Section 
3 we analyze several algorithms for computing the second virial coefficient; this sec- 
tion can be skipped by readers whose main interest is in the results. In Section 4 
we present and analyze our Monte Carlo data for self-avoiding walks in two and 
three dimensions. In Section 5 we compare our results with previous work, discuss 
further the interpretation of the sign of approach to \E'*, and discuss prospects for 
the future. In Appendix A we prove some geometric bounds for subsets of Z ; as a 
corollary, we prove hyperscaling for SAWs in dimension d = 2. In Appendix B we 
discuss the problem of ensuring adequate thermalization for the pivot algorithm. 

®See, for example, the comparisons between theory and experiment in [11, Section 10. F] and [12, 
Section 15.3]. 



7 



In Appendix C we discuss some subtleties involved in the statistical analysis of our 
data. In Appendix D we make a few remarks on the field-theoretic calculations of 
universal amplitude ratios. 



2 Background and Notation 

2.1 The Self-Avoiding Walk (SAW): A Review 

In this section we review briefly the basic facts and conjectures about the SAW 
that will be used in the remainder of the paper. Let C be some regular (/-dimensional 
lattice. Then an N -step self-avoiding walk (SAW) on £ is a sequence of distinct 
points cuo, cui, cujv in C such that each point is a nearest neighbor of its predecessor. 
For simplicity we shall restrict attention to the simple (hyper-)cubic lattice Z ; 
similar ideas apply with minor alterations to other regular lattices. We assume all 
walks to begin at the origin (cuq = 0) unless stated otherwise, and we let Sn be the 
set of all iV-step SAWs starting at the origin (and ending anywhere). 

First we deflne the quantities relating to the number (or "entropy") of SAWs: 
Let Cjv [resp. C]\f{x)] be the number of iV-step SAWs on Z starting at the origin 
and ending anywhere [resp. ending at x]. Then cjv and cn(x) are believed to have 
the asymptotic behavior 

CN ~ /i^iV-^ (2.1) 
cn{x) ~ ^^iV«""8-2 (x flxed ^ 0) (2.2) 

as iV — > oo; here /i is called the connective constant of the lattice, and 7 and Using 
are critical exponents. The critical exponents are believed to be universal among 
lattices of a given dimension d. For rigorous results concerning the asymptotic 
behavior of cjv and cjv(a^), see [1,23,24,98]. 

Next we deflne several measures of the size of an iV-step SAW^: 

• The squared end-to-end distance 

Rl = ul. (2.3) 

• The squared radius of gyration 



(2.4b) 
(2.4c) 




2(iV + 1)^ ,j^o 



^Some other measures of the size of a SAW will be defined in Appendix A. 
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• The mean-square distance of a monomer from the endpoints 



N 

2 



1 

— - y: + - ^nY] . (2.5) 



2(iV + 1) 

We then consider the mean values (i?g)jv, {RI)n and (i?^)jv in the probabihty dis- 
tribution which gives equal weight to each iV-step SAW. Very little has been proven 
rigorously about these mean values, but they are believed to have the asymptotic 
behavior 

{RI)m. {R])n. {RDn ~ N'^ (2.6) 

as iV — > oo, where v is another (universal) critical exponent. Moreover, the ampli- 
tude ratios 

-4« = ^ (2.7) 



N 



Bn = ^ (2.8) 



N 



are expected to approach universal values in the limit N — )• oo.^° Indeed, the full 
probability distribution of cujv is expected to scale as 



N-^"" fixlN"") (2.9) 



as iV — > oo, for a suitable scaling function f (also universal modulo a single nonuni- 
versal scale factor); and / is expected to be rotation-invariant.^^ All these beliefs 
can be subsumed in the even more general assertion that the probability distribu- 
tion of the SAW, with lengths rescaled by iV, converges weakly as iV — )• oo to some 
well-defined probability measure on a space of continuum chains. 

Finally, let cjVi,Af2 be the number of pairs {u^^\u^'^'^) such that cu^^^ is an iVi-step 
SAW starting at the origin, cu^^^ is an iV2-step SAW starting anywhere ^ and cu^^^ and 
cu"^ have at least one point in common (i.e. cu(^) ncu(2) ^ 0). Equivalently, we can 
write cjVi,Af2 in terms of walks that both start at the origin: 

CN,,N. = H T(j'\j'y), (2.10) 



"'^'^ Sometimes the notation = 6(i?^)/(i?g) is usee 



"'^"'^ Actually, (2.9) is claimed to hold only for \x\ of order N'^ . The precise statement of (2.9) is 
therefore that the limit 

exists for each y ^ 0, and that < f(y) < co. 

^^Very recently, Hara and Slade [23,24] have proven that the SAW in dimension d > b converges 
weakly to Brownian motion when N ^ oo with lengths rescaled by CN^I"^ for a suitable (nonuniver- 
sal) constant C . It follows from this that (2.6) holds with = i, and also that (2.7)/(2.8) have the 
limiting values A^o = \, B^o = \- Earlier, Slade [99,100,101] had proven these results for sufficiently 
high dimension d. See also [1]. 



9 



where 

T{J^\J^y) = #{x G I^: J^^ n +x)^0} (2.11) 
is the number of translates of cu^^^ that somewhere intersect co^^K It is beheved that 

CN.,N. ~ ^''^+''HNiN,f^^+^'-'y'g{N,/N,) (2.12) 

as iVi,iV2 — > oo, where A4 is yet another (universal) critical exponent and g is a 
(universal) scaling function. 

The quantity cn^^n^ is closely related to the second virial coefficient. To see this, 
consider a rather general theory in which "molecules" of various types interact. Let 
the molecules of type i have a set Si of "internal states", so that the complete state 
of such a molecule is given by a pair (x, s) where x G Z'^ is its position and s G Si is 
its internal state. Let us assign Boltzmann weights (or "fugacities") Wi(s) [s G Si] 
to the internal states, normalized so that Y^ses, Wi(s) = 1; and let us assign an 
interaction energy Vij (^{x, s), (x' , s' j^ [x,x' G Z'', s G Si, s' G Sj] to a molecule of 
type i at {x,s) interacting with one of type j at (x',s'). Then the second virial 
coefficient between a molecule of type i and one of type j is [102] 



- E E [l-e-^'^^^°''^'^"''''^^] • (2-13) 



s' e s, 



In the SAW case, the types are the different lengths iV, the internal states are the 
conformations io G Sn starting at the origin, the Boltzmann weights are Wjv(^) = 
1/cjv for each io G iSjv, and the interaction energies are hard-core repulsions 

Vnn'[{x,uj),{x ,uj )) = ^ I . ^ ^ (2.14) 

^ ^10 otherwise 



It follows immediately that 



2cjVj CJV2 

The second virial coefficient B^^'^"^"^ is a measure of the "excluded volume" be- 
tween a pair of SAWs. It is useful to define a dimensionless quantity by normalizing 
i?f-^^) by some measure of the "size" of these SAWs. Theorists prefer (Rl) as the 
measure of size, while experimentalists prefer (-R^) since it can be measured by light 
scattering. We thus define the theorists ' interpenetration ratio 

= 2{d/27rY/'^^^ = (d/27rY/'^^^ (2.16) 



and the (usual) interpenetration ratio 

M 



^ 2(rf/12vr)^/^^^ = (rf/12vr)^/^-,^^ (2.17) 
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(for simplicity we consider only Ni = N2 = N). The numerical prefactors are a 
convention that arose historically for reasons not worth explaining here. Crudely 
speaking, \I/ measures the degree of "hardness" of a SAW in its interactions with 
other SAWs. A useful standard of comparison is the hard sphere of radius r and 
constant density: 

R] = -T^r'^ (2.19) 



d + 2 



and hence 



1.12838 md=l 



~ t/r(t//2) ( 3 ) " 1 -1.61859 Zd = l ^^-^^^ 

^2 in (/ = 4 

Inserting (2.1), (2.6) and (2.12) into (2.17), we see that 

^jv ~ iv2^4-7-d- (2.21) 
as iV — > 00. We can therefore distinguish three a priori possibilities: 

(a) 2A4 — 7 — c/i^ > 0, so that \I/jv — > 00 as iV — )• 00. This behavior cannot occur 
unless typical SAWs have a very strange (porcupine-like) shape, which is quite 
implausible.^"^ 

(b) 2A4 — J — diy = 0. In the simplest case this means that \E'jv — > \E'* > 
as iV — )• 00, i.e. typical SAWs exclude each other "within a constant factor 
like hard spheres". This behavior is called hyperscaling . However, the relation 
2A4 — 7 — c/i^ = is also consistent with a logarithmic violation of hyperscaling, 
i.e. \E'jv ~ (log iV)"'' — )• as iV — )• 00 for some power p > 0. 

(c) 2A4 — 7 — c/i^ < 0, so that \E'jv — > as iV — > 00. This is a power-law violation 
of hyperscaling; typical SAWs exclude each other infinitely more weakly than 
hard spheres. 

A very beautiful heuristic argument concerning hyperscaling for SAWs was given 
by des Cloizeaux [103]. Note first from (2.17) that \E' measures, roughly speaking, 
the probability of intersection of two independent SAWs that start a distance of 
order (i?^)^/^ ~ apart. Now, by (2.6), we can interpret a long SAW as an object 
with "fractal dimension" l/v. Two independent such objects will "generically" 
intersect if and only if the sum of their fractal dimensions is at least as large as the 



^^Modulo some reasonable assumptions, this behavior can in fact be rigorously excluded: see 
Theorem A.l and equations (A.20)-(A.22) in Appendix A. 
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dimension of the ambient space. So we expect = lim \I/jv to be nonzero if and 

N^oo 

only if 

- + - > d , i.e. du < 2 . (2.22) 

Since it is believed that 
1 

2 
1 
2 



16 ' 512 

0.588 



for 


(/ > 4 


for 


(/ = 4 


for 


(/ = 4 


for 


(/ = 3 


for 


d = 2 



i + + ^ + for (/ = 4-e [69] (2.23) 



we see 



that 



du <2 for (/ < 4 

c/z/ = "2 + logs" for (/ = 4 
dv > 2 for (/ > 4 



Therefore we expect 



• hyperscaling for d < A 

• logarithmic violation of hyperscaling for d = A 

• power-law violation of hyperscaling for d > A 

One half of this heuristic argument can be proven rigorously. It is easy to see 
[110] that 

CN„N, < (iVl + l)(iV2 + 1)CN,CN, , (2.24) 



so that 



Bf^'""'^ < ^(iVi + l)(iV2 + 1) , (2.25) 



or in terms of critical exponents, 

2A4 - 7 < 2 . (2.26) 

It follows that dv > 2 implies 2A4 — 7 — c/i^ < 0, i.e. the power-law violation of 
hyperscaling. This is now proven rigorously to occur for d > 5 [1,23,24]. 

In the polymer-physics literature it is usually taken for granted that hyperscaling 
holds in dimension d = 3. But in our opinion hyperscaling is a deep property that 
needs to be tested. 

We remark that dimension c/ = 2 is a different case: here hyperscaling can be 
proven rigorously (modulo some reasonable assumptions on the scaling of individual 
SAWs). We present this proof in Appendix A; it is the analogue for SAWs of 
Aizenman's proof [111, Section 8] of hyperscaling for two-dimensional Ising models 
with finite-range ferromagnetic interaction. The underlying geometric idea is that 
SAWs in the plane cannot avoid intersecting each other. 
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Finally, we need to make some comments about corrections to scaling. Clearly, 
(2.1)/(2.2)/(2.6)/(2.12) are only the leading term in a large- iV asymptotic expan- 
sion. According to renormalization-group theory [112], the mean value of any global 
observable O behaves as iV — )• oo as 



{0)n = AN^ 



1 + — H - + ...-{ — H \ = h . . . 

, Co Ci C2 

' 7\TA_I1 ' 7\TA_IO~r*** 



(2.27) 



Thus, in addition to "analytic" corrections to scaling of the form ak/N^ , there are 
"non- analytic" corrections to scaling of the form bk/N^^~^^ , Ck/N^^~^^ and so forth, 
as well as more complicated terms [not shown in (2.27)] which have the general form 
const / N^^^^~^^^^^~^ where ki,k2, ■ ■ ■ and / are non-negative integers. The leading 
exponent p and the correction-to-scaling exponents Ai < A2 < . . . are universal; p 
of course depends on the observable in question, but the do not. [Please note 
that the exponents Ai < A2 < . . . have no relation whatsoever to the gap exponent 
A4 defined in (2.12). The notation used here is standard but unfortunate.] The 
various amplitudes (both leading and subleading) are all nonuniversal. However, 
ratios of the corresponding amplitudes A, 60 and Cq (but not or the higher b^^ Ck) 
for different observables are universal [88,13]. 

Remark. The names of the critical exponents 7, asm^, ^ and A4 are chosen by 
analogy with the corresponding exponents in ferromagnetic spin systems [17,113]. 
Indeed, the generating functions of self-avoiding walks, 

00 

x(/3) = E/^'^^iv (2.28) 

iV=0 

00 

G{x-P) = E/3^c4x) (2.29) 

N=0 

00 

U4(/3) = -3 E 13'''^''' CM„N, (2.30) 

Ni,N2=0 

are equal to the susceptibility, spin-spin correlation function and fourth cumulant 
in the n-vector model analytically continued to n = [1,4,5,6,7]. In particular, if x 
is a nearest neighbor of the origin, then G(x; fi) is essentially the energy E (up to 
an additive and multiplicative constant). The quantity 




m ^ (2.31) 



is the second-moment correlation length. Inserting (2.1)/(2.2)/(2.6)/(2.12) into 
(2.28)-(2.31), we obtain the leading behavior 

x(/3) ~ iPc-^r' (2.32) 
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m 

M4(/3) 



(/3c - ^r^ 



+ regular terms 



(2.33) 
(2.34) 
(2.35) 



as /3 approaches the critical point fic = Note, in particular, that a^ing is the 

exponent for the singular part of the specific heat Ch ~ dE/djS; the exponent for 
the full specific heat is a = niax(asing^ 0). If the hyperscaling relation = 2A4 — 7 
holds (without multiplicative logarithmic corrections), then the renormalized cou- 
pling constant g = —U4^/x^^'^ tends to a nonzero limiting value </* as /3 — )• [ic, so 
hyperscaling (without multiplicative logarithmic corrections) can be interpreted as 
the non-Gaussianness (nontriviality) of the scaling-limit quantum field theory [7, 
p. 281]. 

2.2 The Pivot Algorithm: A Review 

The pivot algorithm was invented in 1969 by Lai [8], reinvented in 1985 by 
MacDonald et al. [9], and again reinvented a short time later by Madras [10]. The 
pivot algorithm is the most efficient algorithm currently known for simulating SAWs 
in the fixed- iV, variable-x ensemble. Here we summarize briefly the relevant features 
of the algorithm; more details can be found in [10,81]. 

The elementary move of the pivot algorithm is as follows: Choose at random a 
pivot point k along the walk (0 < k < N — 1); choose at random an element g of the 
symmetry group of the lattice (rotation or reflection or a combination thereof); then 
apply g to the part of the walk subsequent to the pivot point (namely iOk+i^ • • • , ^jv), 
using iOk as the temporary "origin". That is, the proposed new walk io' is 



The walk io' is accepted if it is self-avoiding; otherwise, it is rejected, and the old walk 
id is counted once more in the sample. It is easy to see that this algorithm satisfles 
detailed balance for the standard equal-weight SAW distribution. Ergodicity is less 
obvious, but it can be proven [10,114]. 

At flrst thought this seems to be a terrible algorithm: for N large, nearly all the 
proposed moves will get rejected. In fact, this latter statement is true, but the hasty 
conclusion drawn from it is radically false! The acceptance fraction / does indeed 
go to zero as iV — )• 00, roughly like iV"''; empirically, it is found that the exponent p 
is sa0.19 in (/ = 2 [10] and saO.ll in d = 3 [10,115,116]. But this means that roughly 
once every moves one gets an acceptance. And the pivot moves are very radical: 
one might surmise that after very few accepted moves (say, 5 or 10) the SAW will 
have reached an "essentially new" conflguration. One conjectures, therefore, that 
the autocorrelation time r of the pivot algorithm behaves as Things are 

in fact somewhat more subtle (see the next paragraph), but roughly speaking (and 
modulo a possible logarithm) this conjecture appears to be true. On the other hand, 
a careful analysis of the computational complexity of the pivot algorithm (see also 
below) shows that one accepted move can be produced in a computer time of order 




for < z < A; 

for A; + 1 < z < iV 



(2.36) 
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N. Combining these two facts, we conclude that one "effectively independent" 
sample can be produced in a computer time of order N (or perhaps NlogN). 

Let's look more closely: Suppose we know that the acceptance fraction / in the 
pivot algorithm behaves as / ~ iV"'' as iV — )• oo. Then, as argued above, after a 
few successful pivots — i.e. a time of order 1// ~ — the global conformation of 
the walk should have reached an "essentially new" state. Thus, we expect that for 
observables A which measure the global properties of the walk — such as or 
i?^ — the autocorrelation time Tint,A [see (C.4)] should be a few times 1//. This 
is confirmed numerically [10, Section 4.3]. On the other hand, it is important to 
recognize that local observables — such as the angle between the 17^^ and IS*'' steps 
of the walk — may evolve a factor of N more slowly than global observables. For 
example, the observable mentioned in the preceding sentence changes only when coir 
serves as a successful pivot point; and this happens, on average, only once every N/ f 
attempted moves. Thus, for local observables A we expect Tint,A to be of order N/ f. 
General properties of reversible Markov chains [77] then imply that the exponential 
autocorrelation time Texp must be of at least this order; and if we have not overlooked 
any slow modes in the system, then Texp should be of exactly this order. Finally, 
even the global observables are unlikely to be precisely orthogonal to the slowest 
mode; so it is reasonable to expect that Texp,A be of order N/ f for these observables 
too. In other words, for global observables A we expect the autocorrelation function 
PAA{t) to have an extremely-slowly-decaying tail which, however, contributes little 
to the area under the curve. This behavior is illustrated by the exact solution of 
the pivot dynamics for the case of ordinary random walk [10, Section 3.3], and by 
numerical calculations for the SAW (see Appendix C). 

Computational complexity. A very important issue in any algorithm — but 
especially in a non-local one — is the CPU time per iteration. By using a hash 
table [117,118], the self-avoidance of a proposed new walk can be checked in a time 
of order N . But one can do even better: by starting the checking at the pivot point 
and working outwards, failures can be detected in a mean time of order iV^"'' [10, 
Sections 3.4 and 4.4]. The mean CPU time per successful pivot is therefore ~iV^~'' 
for each of ~iV failures, plus ~iV for one success, or ~iV in all. Combining this with 
the observations made previously, we conclude that one "effectively independent" 
sample — as regards global observables — can be produced in a computer time of 
order N . 

Initialization. There are two main approaches: 

1) Equilibrium start. Generate the initial configuration by dimerization [119,1,81]; 
then the Markov chain is in equilibrium from the beginning, and no data need be 
discarded. This approach is feasible (and recommended) at least up to N of order 
a few thousand. There is no harm in spending even days of CPU time on this 
step, provided that this time is small compared to the rest of the run; after all, the 
algorithm need only be initialized once. 

2) "Thermalization". Start in an arbitrary initial configuration, and then dis- 
card the first Udisc ^ Texp ~ N/f iterations. This is painful, because Texp is a factor 
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~iV larger than Tint,A for global observables A; thus, for very large N (> 10^), the 
CPU time of the algorithm could end up being dominated by the thermalization. 
Nevertheless, one must resist the temptation to cut corners here, as even a small 
initialization bias can lead to systematically erroneous results, especially if the sta- 
tistical error is small: see Appendix B for striking evidence of this. Some modest 
gain can probably be obtained by using closer-to-equilibrium initial configurations 
(e.g. [120]), but it is still prudent to take ridisc at least several times N/f. 

Initialization will become a more important issue in the future, as faster com- 
puters permit simulations at ever-larger chain lengths. 

3 Algorithms for Counting Overlaps 

In this section we discuss algorithms for computing the excluded volume between 
a given pair of SAWs; this is the key step in a Monte Carlo study of the second 
virial coefficient. This section can be skipped by readers whose main interest is in 
the results rather than the algorithms. 

3.1 Generalities 

Let io^^'^ and cu^^^ be, respectively, iVi-step and iV2-step SAWs, and define T(c<;^"^\ cu^^^) 
to be the number of translates of cu^^^ which somewhere intersect cu^^^: 

T{J^\J^^) = #{x G Z'^: cu(^) n + x) 7^ 0} (3.1a) 
= #(cu(i) - cu(2)) , (3.1b) 

where A — B = {y — z: y A^ z B}. The expected value of T(c<;(^\ cu^^^), averaging 
over independent walks co^^'^ G Sn^ and cu^^^ G Sn^^ is cn^^n^/cn^cn^. This quantity 
has the asymptotic behavior 

cnunJcn.cn, ~ imN2 f^^-'^^/'g{m/N2) , (3.2) 

where g is a scaling function [cf. (2.1)/(2.12)]. It is thus possible to estimate the 
critical exponent 2A4 — 7 by running two independent pivot algorithms and mea- 
suring T((jj^^\(jj^'^'^). [Typically one would run at Ni = N2 = N for a sequence of 
values of N .] In particular, this allows a direct Monte Carlo test of the hyperscaling 
relation = 2A4 — 7. Note that an independent measurement of 7 is not needed. 

The efficient determination of T(Ai,A2) = i^{Ai — A2) for a specified pair of 
sets Ai, A2 C Z is a very interesting and nontrivial problem in computer science. 
We see two broad approaches: 

1) Deterministic algorithms which compute T(Ai, A2) exactly. 

2) Monte Carlo algorithms which produce an unbiased (or almost unbiased) es- 
timate of T(Ai, A2). 
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In the latter case, the statistical fluctuations in the auxiliary (inner-loop) Monte 
Carlo process would be added to those in the main Monte Carlo program; but this 
is acceptable provided that the former are not too large compared to the latter (see 
Section 3.4). 

We shall discuss the deterministic algorithms for computing T(Ai, A2) in Section 
3.3, and the Monte Carlo algorithms in Sections 3.4-3.9. 



3.2 Notation 

We shall denote by J\fi (resp. A/2) the number of points in the set Ai (resp. 
A2). We also write Afmin = min(A/i,A/2) and Afmax = max(A/i, A/2). An iV-step 
self-avoiding walk has J\f = N -\- 1 points. 

Now flx a pair of sets Ai, A2 C Z'', and write 

S = Ai- A2 = {y - z: y e Ai, z e A2} . (3.3) 

Our goal is to compute T(Ai,A2) = i^{S). An important role is played by the 
function 

Pix) = #{{y,zy. y e Ai, z e A2, y - z = x} (3.4a) 
= iXA^ * XA2 )(x) , (3.4b) 

where xAi and XA2 ^.re the indicator functions of the sets Ai and A2, respectively, 
and XA2{^) = XA2{~^)7 ^.nd * denotes convolution. Clearly 

p{x) = for X (ji S (3.5a) 

1 < p{x) < Mmin for X ^ S (3.5b) 



We also write 



Note that 



Y^Hx) = El = T{A,.A2) (3.7) 
E P{^) = E Pi^) = MAr2 (3.8) 



For future reference we deflne also 



C.(-4„A,, .Z-i^. ,3.9) 
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This observable has httle intrinsic interest, but it wiU play an important role in 
two of the Monte Carlo algorithms. It is not hard to see that^"* 

A/'i+A/'2-l < T(Ai,A2) < A/'iA/'2 (3.10) 

^^^^^j^f < U{AuA,) < T{A,,A,) (3.11) 
Moreover, we shall prove in Appendix A (Theorem A. 3) that 

UiA^,A2) > log(MAr2) . (3.12) 

We shall use ( • ) to denote expectation with respect to some probability distri- 
bution on pairs of sets Ai, A2 (e.g. the equal- weight ensemble on the space Sn^ X Sn^ 
of pairs of SAWs). We shall use E( ■ ) to denote expectation with respect to some 
"inner-loop" Monte Carlo algorithm. 



3.3 Deterministic Algorithms 

Here we introduce some deterministic algorithms for computing p{x)^ 
T{Ai, A2) and/or U{Ai,A2). These algorithms can be employed either "stand- 
alone" or as building blocks for the Monte Carlo algorithms to be introduced later. 



1) Suppose first that we want to compute p{x) for a single value of x. This 
can be done as follows: write all the points of Ai into a hash table [118,117]; then 
examine sequentially each of the points 2; G A2, inquiring whether y = x-\- z belongs 
to Ai, and incrementing a counter if it does. Clearly this requires a CPU time of 
order J\f\ -\- A/2 . 

For computing /(x), the algorithm can be streamlined by stopping as soon as 
one finds p{x) > 0. 

2) Suppose next that we want to compute p{x) for several different values of x, 
say r of them. Then the foregoing algorithm requires a CPU time of order J\fi -l-rA/2. 
But since the answer is invariant under the interchange Ai A2, x — )• — x, it now 
pays to choose A2 to be the smaller of the two sets. As a result, the CPU time is 
of order Af^ax + rAf^^n- 

3) Suppose, finally, that we want to compute p{x) or I{x) for all x. [A special 
case of this is to compute T(Ai, A2) or U{Ai, ^2)-] This can be done by examining 
sequentially each of the pairs y G Ai, 2; G A2, and writing the points x = y — z into 
a hash table equipped with an auxiliary count field. Clearly this requires a CPU 
time of order A/1A/2. [The CPU time for subsequently utilizing the counts {p{x)} is 

^■*The two upper bounds are trivial. The lower bound in (3.11) is the Schwarz inequality. The 
lower bound in (3.10) is proven as follows: Let j/* (resp. z*) be the lexicographically smallest element 
of Ai (resp. A2). Then y^, — A2 ^ y^, — z^, ^ Ai — z^,; so j/* — A2 and Ai — are subsets of S (of 
cardinalities Af2 and Afi, respectively) with only one point in common (namely, j/* — 2*)- Q.E.D. 
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of order T(Ai, A2) < A/1A/2.] If all one wants is {/(x)}, then the count field can be 
dispensed with; it suffices to know which sites x are hit at least once. 

4) An alternative algorithm for computing p{x) for all x can be based on the 
Fast Fourier Transform (FFT). Let aj (1 < j < d) be the extension of Ai in the j^^ 
coordinate direction (i.e. the difference between the maximum and minimum values 
of j/j for y G Ai). Let /3j be the corresponding extension for A2. Now let kj be the 
least integer such that 2^^ > aj + /3j ■ + 1- Then if we place Ai and A2 in a periodic 
box B of size 2^^ X 2*^^ x • • • x 2*^'*, we can compute the convolution of xAi and XA2 
in the box 5, without distortion by the periodic boundary conditions; this gives 
the full set of counts {p{x)} [and thus also {/(x)}, T(Ai, A2) and U(Ai^ A2)]. The 
convolution can be carried out by the FFT in a time ~ V logV ^ where 

d 

V = 2'''+-+'''' < l[2(aj+/3j). (3.13) 

Let us now estimate the performance of these methods as "stand-alone" algo- 
rithms for computing T(c<;(^\ cu^^^), where cu^^^ and cu^^^ are self-avoiding walks. Al- 
gorithm #3 computes T(cu(^), cu^^)) in a CPU time of order MA/'2 = (Ni + 1)(N2 + 1), 
i.e. of order iV^ if Ni ~ N2 ~ N . By contrast, we expect that one "effectively inde- 
pendent" sample of the pair (cu^^^ cu^^^) can be produced by the pivot algorithm in 
a CPU time of order Ni -\- N2 (if Tint,T ~ N^) or in any case not much greater. So 
this algorithm would spend more time analyzing the data than producing it! — and 
the overall computational complexity per "effectively independent" sample would 
be increased from N to iV^, thereby nullifying the advantage of the pivot algorithm 
over previous [121,122] algorithms. 

Algorithm ^4 computes T((jJ^^\(jJ^'^'^) in a CPU time of order UlogU, where V 
is given by (3.13). For "typical" SAWs we have V ~ (iVi -|- N2Y'^ , so presumably 
we have 

{V log V) ~ (iVi + N2Y'' log(iVi + N2) . (3.14) 

In the usual situation Ni ~ N2 ~ iV, this method is asymptotically better than 
algorithm ^3 if < 2, i.e. if d < 4. However, it is unlikely to be better in practice 
except for very large N. And the behavior is still vastly worse than the time of 
order N for generating the walks (except in d = 1). 

It may be possible to devise deterministic algorithms which are more efficient 
than either of these elementary ones; we leave this as an exercise for interested 
computer scientists. 

3.4 Monte Carlo Algorithms: Generalities 

Fix a pair of sets Ai, A2, and suppose that we use some Monte Carlo algorithm 
to provide an unbiased estimate Z of T(Ai, A2). We will thus have 

E{Z\Ai,A2) = T(Ai,A2) (3.15a) 
var(Z|Ai,A2) = E{Z^\A^, A2) - E{Z\A^, A2Y = ^(^1,^2) (3.15b) 
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Here (3.15a) expresses the unbiasedness of the inner-loop Monte Carlo algorithm, 
while (3.15b) defines its (conditional) variance. We shall compute the functional 
y(Ai, A2) for each of the Monte Carlo algorithms we introduce (Sections 3.5-3.7). 

Now let us call the Monte Carlo subroutine R times for the given pair (^1,^2) 
and average the results: Z = YliLi ^i- Obviously we have 

E{Z\Ai,A2) = T(Ai,A2) (3.16a) 
var(Z|Ai,A2) = R-^V{Ai,A2) (3.16b) 

Now suppose that we generate a random pair (^1,^2) from some probability 
distribution [e.g. SAWs ((jJ^^\ io^'^'^) from the equal-weight distribution on Sn^ X Sn^]. 
Clearly Z is an unbiased estimator of (T), i.e. 

{Z) = {EiZ\A„A2)) = (T) , (3.17) 

where ( • ) denotes expectation in the given probability distribution. The variance 
of Z is a sum of two terms: 

var(Z) = var(T) + (var(Z|Ai, A2)) 

= [(T^) - (T)2] + R-^V) . (3.18) 

The first term is the fiuctuation of T(Ai, A2) from one pair of sets to another; the 
second term is the mean over pairs of sets of the fiuctuation (conditional variance) 
in the inner Monte Carlo subroutine. 

The mean CPU time for the computation of Z is (Tcpu) = a + bR: here a is 
the mean CPU time for generating a pair of "effectively independent" sets (Ai, A2) 
from the desired ensemble, plus any "setup" time associated with the inner-loop 
Monte Carlo algorithm; while b is the mean additional CPU time per iteration of 
the inner-loop Monte Carlo algorithm. The goal is to minimize the variance-time 
product 



{Tc 



^cpu)yaT(Z) = bvaT(T)R + a{V)R-^ + avar(T) + b{V) , (3.19) 

since this quantity divided by the total CPU time equals the variance of our final 
estimate. Hence the optimal choice of R is 

and the variance-time product is then 

[{TcPu)y^T{Z)U = [(avar(T))i/2 + {b{Vyf'f . (3.21) 



Of course, R must be a positive integer, and so the true Ropt is obtained by rounding 
the right-hand side of (3.20) up or down. This subtlety can be ignored if Ropt is 
large, but may be significant otherwise. In particular, the deterministic inner-loop 
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algorithms (Section 3.3) have V = 0^ but in this case Ropt = 1 [rather than as 
(3.20) claims] and [{Tcpu) ^s.T(Z)]opt = (a + 6)var(T). 

In the remainder of this section we will assume that the CPU time for generating 
the sets Ai and A2 is of order J\fi + A/2. Clearly this is a lower bound, since it takes a 
time of order J\fi + A/2 simply to write down the two sets. On the other hand, for our 
application to SAWs, Ai and A2 will be generated by the pivot algorithm, which 
generates an "effectively independent" SAW (as regards global observables) in a 
CPU time of order N (see Section 2.2). We expect that the observable T((jJ^^\(jJ^'^'^) 
is indeed "global" in the sense that Tint,T ~ where p is the acceptance-fraction 
exponent. 

3.5 "Hit-or-Miss" Monte Carlo Algorithm 

Let aj' (resp. aj) be the maximum (resp. minimum) value of the j^^ coordinate 
among the points in Ai, so that 

B, = [a-, at] X ... X (3.22) 

is the smallest rectangular parallelopiped containing Ai. Let [3^ and /3~ be the 
corresponding values for A2, and B2 the corresponding box. It follows that 

B = B^-B2 = [a- - a+ - X . . . X [a," - a+ - (3.23) 

is a parallelopiped which is guaranteed to contain all the points of S" = Ai — A2^^ 
Therefore, T(Ai,A2) = i^{S) = Y^xeB ^{^) can be computed by the trivial "hit 
or miss" Monte Carlo method: Pick a point x ^ B at random, and compute I{x) by 
the deterministic algorithm ^1 of Section 3.3; then output Z = jj^{B) I{x), where 
#(5) = 11^=1 ("j + 1^3 + !)• [Here = - aJ and /3j = /3+ - /3~ .] Clearly I{x) is 
a binomial random variable of mean p = ^(S)/^(B)^ so that 

Eim) = Illy (3.24) 

var(/(x)) = (3.25) 

Hence Z = jj^{B) I{x) is an unbiased estimator of ^(S"), and its variance is 

VHH-or-^^ss{Ar,A2) = var(Z) = #(5) [#(5) - #(5)] . (3.26) 

The CPU time for R iterations of this algorithm is of order Tcpu = J^max + 
RJ^min'- we put the larger of the two sets Ai, A2 in the hash table once, and then 
each time we compute p{x) by looping over the smaller of the two sets. The CPU 
time for generating the two sets is by assumption also of order Mmax- Therefore, in 
(3.19)-(3.21) we have a ~ Mmax and h ~ Mmin- 



^®In fact it is the smallest such parallelopiped, since for each index j there is a point in Ai — A2 
with j*^ coordinate equal to a~ — (3f , and another point with j*'' coordinate equal to — (3~ ■ 
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3.6 Barrett Algorithm: Theory 

Barrett [123] has proposed the foUowing Monte Carlo algorithm, which gives an 
unbiased estimate of T(Ai, A2): 

1. Choose at random y G Ai and z ^ A2. Set x = y — z. 

2. Compute p{x) using the deterministic algorithm ^1 described in Section 3.3. 
[Note that by construction we have x ^ S and hence p{x) > 0.] 

3. Output Y = MiM2lp{x). 

The analysis of this algorithm is easy: In Step 1 we choose the vector x with 
probability 

p{x) 



It follows that 



E{Y) = ^Prob(x)F(x) = ^1 = T(Ai,A2) (3.28a) 

E(Y') = ^Prob(x)F(x)2 = = ^i^2U(A,,A2) (3.28b) 

xes xes P\^) 

and hence 

VBarrettiA^,A2) = Var(F) = Af^Af2U (A^ , A2) " ^^,^2)^ . (3.29) 

The CPU time for R iterations of the Barrett algorithm is of order Tcpu ~ 
J^max + RJ^min'- we put the larger of the two sets Ai, A2 in the hash table once, 
and then each time we compute p{x) by looping over the smaller of the two sets. 
The CPU time for generating the two sets is by assumption also of order Mmax- 
Therefore, in (3.19)-(3.21) we have a ~ Mmax and h ~ Mmin- 

3.7 Karp-Luby Algorithm: Theory 

Karp and Luby [124,125] have devised an elegant Monte Carlo algorithm for 
estimating T{Ai, A2) [and somewhat more general combinatorial problems]. The 
Karp-Luby algorithm goes as follows: 

1. Choose at random y G Ai and z ^ A2. Set x = y — z. Set t = 1. 

2. Choose at random y' G Ai. 

3. If 2;' = y' — x G A2, then go to Step 4. Otherwise, increment t by 1 and go to 
Step 2. 

4. Output Z = N'2t. 
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This algorithm can be understood as a randomized version of the Barrett algorithm. 
Step 1 is identical in the two algorithms, and it selects the vector x with probability 

Prob(x) = ■ (3.30) 

Then t (the number of trials of Steps 2 and 3 needed to find a y' such that y' — x G A2) 
is, conditioned on x, a random variable with a geometric distribution: 

Prob(t = k\x) = ^ (^1 - for k = 1,2,3,... . (3.31) 

Hence the conditional expectations of Z = M^^^ are 

E{Z\x) = (3.32a) 
p(x) 

^ ^ - 1) (3.32b) 

(Of course, this makes sense only for x G S.) Thus, Steps 2-4 of the Karp-Luby 
algorithm produce a random quantity Z whose mean value (conditional on x) is 
precisely the deterministic quantity Y = J\fiJ\f2 / p{x) of the Barrett algorithm. It 
follows that the unconditional expectations are 

E{Z) = J2FToh{x)E{Z\x) = T(Ai,A2) (3.33a) 

xes 



E{Z') = ^Prob(x)E(Z2|x) = AT^^ [^-1' 



x^S x^z 

2A/'iA/'2t/(Ai,A2) - A/'2T(Ai,A2) 

(3.33b) 



and hence 



VKarp-LubyiAi,A2) = Var(Z) = 2AfiAf2U (A^ , A2) - Af2T{Ai,A2) - T{Ai,A2f . 

(3.34) 

For future reference, we note the following inequality: 

' \ 1 ) -^l-^2^( Ai , A2) < VKarv-Luhy{Ar,A2) < 2M^M2U (A^ , A2) . 

(3.35) 

[Proof: The upper bound is trivial. To prove the lower bound, use (3.10) and 
(3.11) to deduce (A/i + A/2 — 1)T < < A/iA/2?7.] This means that VKarp-Luby is 
of the same order of magnitude as its first term, namely ~ A/iA/2?7, except perhaps 
when J\fi <^ A/2. An alternative (and often sharper) lower bound on VKarp-Luby can 
be obtained by noting that 

^2T < (^^logAC.^.) \^A'2U . (3.36) 
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[Proof: From (3.11) we have T'^ / M1M2 < U, while from (3.12) we have logJ\fmax < 
U . Now take the geometric mean of these two bounds.] Therefore we have 



VKa 



rp — Luby 



(Ai,A2) > 



log A/'ma. 



-1/2 



MAr2t/(Ai,A2) - T{A^,A2f 



(3.37) 

The CPU time for one execution of Steps 2 and 3 is essentially t, so the expected 
CPU time for one iteration of the Karp-Luby algorithm is E(t) = T(Ai, A2)/A/2. In 
addition, there is an initial CPU time of order A/2 to place the elements of A2 in a 
hash table. Finally, we should remember the time of order J\fi + A/2 for generating Ai 
and A2 in the first place. The expected CPU time for R iterations of the Karp-Luby 
algorithm (plus generating Ai and A2) is thus 



r(Ai,A2) 

I CPU ~ J^max H 77 -K 

J\2 



(3.38) 



Therefore, in (3.19)-(3.21) we have a ~ Afmax and b ~ (T)/A/2. 



3.8 Scaling Theory 

In this section we consider the scaling theory of the three Monte Carlo algorithms 
— hit-or-miss, Barrett and Karp-Luby — in the case where Ai and A2 are indepen- 
dent random SAWs of lengths Ni and iV2, respectively, and Ni ~ N2 ~ iV — )• 00. In 
each case we need to compute (or guess heuristically) the scaling behavior of var(Z) 
and (Tcpu) as iV — )• 00. A good figure of (de)merit for an algorithm is the mean 
CPU time needed to estimate (T) with a relative variance of order 1; this time is 



[{Tcpu) var(Z)]opt 



(3.39) 



For the hit-or-miss algorithm, we must study the scaling of (T), (T^) and (T ^(B)). 
For the Barrett and Karp-Luby algorithms, we must study the scaling of (T), (T^) 
and (U) as well as their various combinations. 

As discussed in Section 2.2, we expect that (T) scales as 

(T) ~ , (3.40) 

where 

PT = 2A4 - 7 . (3.41) 

We further expect that the probability distribution of T will, after rescaling by 
iV^, approach a nontrivial limiting distribution (here "nontrivial" means that the 
distribution is not a delta function). Therefore, we expect that 

(T^), var(T) ~ iV^f^ . (3.42) 

On the other hand, it is reasonable to expect that 

{U) ~ N'P^ (3.43) 
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for some (a priori unknown) exponent pu- From (3.11) we know that 2px — 2 < 
Pu 1^ Pt] and we will be interested in knowing whether these inequalities are strict 
or not. 

We now look at the individual algorithms: 

Hit-or-Miss Algorithm: As just discussed, we expect 

(T) ~ (3.44) 
var(T) ~ N"^"^ (3.45) 

with pt = 2A4 — 7. On the other hand, for "typical" pairs of SAWs we clearly have 
#(5) ~ N"^", so we expect 

{T#{B)) ~ N'^'^+'P^ . (3.46) 

In general we have > px] but even when = px (i.e. hyperscaling holds), it 
seems intuitively clear that the ratio (T^)/(T ^(5)) will stay well below 1 (except 
in dimension d = 1). Therefore, radical cancellations in (3.26) are excluded, and we 
expect 

{VMt-or-m^ss) ~ N'^+P^ . (3.47) 

It follows from (3.45) and (3.47) that 

var(Z) ~ N"^"^ + R-^N'^^+p^ . (3.48) 
On the other hand, the CPU time is obviously 

{Tcpu) ~ N + RN . (3.49) 
The variance-time product is thus 

(Tcpc/) var(Z) ~ n[n'^p^R + N'^^+p^ R-^ + n'^pt + p^dr^+prj _ ^g^^Q^ 

There are two cases: 

(a) If the hyperscaling relation dv = px holds, then Ropt ~ 1. [We expect that 
this is the case for d < 4.] With this choice of i?, we have var(Z) ~ (T)^ and 
(Tcpu) ~ N . In other words, we can achieve a relative variance of order 1 in a CPU 
time of order N . 

(b) If du > Pt, then Ropt ~ N(d'^-PT)/2 _ expect that this is the case for 
d > 4. More precisely, for d > 4 we expect that = 1/2 and px = 2, in which 
case Ropt ~ N('^-*y*.] With this choice of R, we have var(Z) ~ (TYI'^ N'^"!'^ and 
(Tcpu) ~ iV^+''''/^(T)~^/^, which implies a variance-time product of order N^'^'^'^ (T) . 
[In fact, a variance-time product of this order can be achieved with any R in the 
range 1 < i? < N'^'^~'''^ ; the optimal value Ropt lies at the geometric mean of these 
two extremes.] Hence a relative variance of order 1 requires a CPU time of order 
p^i+du-PT [i_g_ order N^+(d-4)/2 for d > 4]. 
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Barrett and Karp-Luby algorithms: The behavior of the Barrett and Karp- 
Luby algorithms is determined by the scahng behavior of (T), (T^) and {U) as 
N ^ oo. From (3.29) and (3.35)-(3.37), we have 



(VBarrett) = M^M2{U) 



{'^Karp — Luby 



1 



MrM2{U) 
MiM2{U) 



(3.51) 



(3.52) 



[We write a 6 to denote that a/6— )-lasiV— )-oo.] This focusses attention on the 
ratio 

[which by (3.11) satisfies < n{Ni,N2) < 1] and in particular on 

K = lim k{N,N) . (3.54) 

N^oo 

(For simplicity we consider here only Ni = N2 = iV, but the same principles ob- 
viously apply at any fixed ratio N1/N2 7^ 0,oo.) There are then three possible 
cases: 

(a) K = 0. In this case {VBarrett) ~ N^{U) and {VKarp-Luby) ~ 2N^{U). 



(b) < /€ < 1. In this case (VBarrett) ^ (1 - k)N'^{U) and {VKarp-Luby) ^ (2 - 

k)N^(U). 

(c) n = 1. In this case {VBarrett) ^ N'^{U), and its exact scaling is subtle. Of 
course we still have {VKarp-Luby) ~ N'^{U). 

Case (a) corresponds to the exponent pu being strictly greater than 2px — 2, while 
cases (b) and (c) correspond to pu = 2pT — 2. We believe that in fact case (c) never 
occurs: although it is possible to have U = T^/A/iA/2 for special pairs (jJ^^\ (jJ^'^\ e.g. 
perpendicular rods, it seems quite implausible that such behavior could occur (even 
in the limit N — )• 00) after averaging over all pairs of SAWs. 

In Section 3.9 we shall present numerical data showing clearly that «: < 1 in 
dimensions c/ = 2, 3; what is less clear is whether k, is zero or nonzero. For dimension 
d > 4 we expect that < k, < 1. [For d > 4 we expect that (T) ~ iV^/^ from which 
it follows that {U) ~ and hence < k < 1}^ On the other hand, we believe, as 
noted above, that «: = 1 is impossible.] 

Assuming that < /€ < 1, it follows that {VBarrett) and {VKarp-Luby) are both of 
order N'^{U). Therefore, 

var(Z) ~ (T)2 + R-^N^U) (3.55) 



^^This is rigorously proven for > 6, modulo the problem of translating generating-function 
results into fixed-# results. See [1, Theorem 1.5.5 and Remark following it]. 

i^Proof: {U) > {T^)/iN + if > {Tfl{N + 1)^ by (3.11) and the Schwarz inequality; while 
{U) < (T) by (3.11). Thus, if (T) ~ we have also {U) ~ N"^ . 
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for both algorithms. On the other hand, the CPU times are 

(TcPU,Barrett) ~ N + NR (3.56a) 
{TcPU,Karp-Luby) ~ N + N-^T)R (3.56b) 

The variance-time products are thus 

{TcPU,Barrett)^^r{Z) ~ N\{TYR + {U) R-' + (T)^ + {U)] 

(3.57a) 

{TcPu,Karp-Luby)y^r{Z) ~ n[n-'{T)^R + N'{U)R-' + (T)^ + (t/)(T)] 

(3.57b) 

RecaUing now the inequahty < (N + 1)^?7 [cf. (3.11)], we see that the dom- 
inant term for the Barrett algorithm is the fourth one, provided that 1 < i? < 
N'^ (U) / {Ty ; and in this case we have (Tcp[/ ) var(Z) ~ N^{U). The optimal value 
is at the geometric mean of this range, i.e. Ropt ~ N{Uy/y{T). We conclude that 
a relative variance of order 1 requires a CPU time of order (U) / {T^ . This is 
at least of order N; but it may be larger, in case (U) ^ (T^)/iV^ (i.e. in case 
pu > 2pT - 2). 

Recalling next the inequality U <T [cf. (3.11)], we see that the dominant term 
for the Karp-Luby algorithm is the third one, provided that N'^ (U) / {T^ ^ R ^ 
N'^/{T); and in this case we have ( Tcp[/) var(Z) ~ iV(T)^. The optimal value is at 
the geometric mean of this range, i.e. Ropt ~ iV^(?7)^/^/(T)"^/^. We conclude that a 
relative variance of order 1 requires a CPU time of order N . 

Comparing the analyses of the Barrett and Karp-Luby algorithms, we see that 
the effect of the additional randomization in the Karp-Luby algorithm is to reduce 
drastically the mean CPU time per iteration (T/A/2 versus A/'mm) while only mod- 
estly increasing the variance. The Karp-Luby algorithm is thus superior to the 
Barrett algorithm whenever (T) <^ N"^ (as we expect occurs in dimension c/ < 4); 
the two algorithms are of the same order whenever (?7) ~ (T) ~ iV^ (as we expect 
occurs in dimension d > A). 

3.9 Barrett and Karp-Luby Algorithms: Numerical Results 

In this section we report results of a rather crude Monte Carlo study of (T), 
(T^), var(T), (?7), {VBarrett) and {VKarp-Luby) for SAWs in dimensions d = 2 and 
d = 3 (taking Ni = N2 = N). The goal is to estimate the exponents pp and pu 
(and the various amplitudes) in the scaling relations 
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r-^ 


ApNP^ 


(3.58a) 








(3.58b) 


var(T) 






(3.58c) 
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(3.58d) 
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In Table 1 we report our Monte Carlo estimates for 2-dimensional self-avoiding 
walks at iV = 100, 200, 300, 400, 500. These were obtained by using the dimer- 
ization algorithm to generate pairs of independent SAWs cu^^^cu^^^ and then using 
deterministic algorithm ^3 (Section 3.3) to compute T((jJ^^\(jJ^'^'^) and U((jJ^^\(jJ^'^'^). 
Unfortunately, some of the observables lack error bars (this is our fault). The fit 
to (T) gives a reasonable provided that we discard the data point at the lowest 
value of N; we then get px = 1.494 it 0.001. This is not far from the expected 
value pt = = 3/2; the small discrepancy can quite plausibly be attributed to 
corrections to scaling. (We will make a much more careful estimate of (T) and its 
exponent px in Section 4.1.) The fits to (T^) and var(T) are also consistent with 
Pt = 3/2, although the lack of error bars prevents making this quantitative. The 
fit to (VKarp-Luby) gives a reasonable only if we discard the data points at the 
two lowest values of N; we then get pu = 1.110 it 0.002. If we use the same data 
points in the fits to (U) and (Vearrett), we get pu = 1.078 and 1.152 (unfortunately 
without error bars). These estimates of pu are thus in rather mediocre agreement; 
and it is far from clear whether pu is strictly greater than 2px — 2 = 1. Another way 
of looking at this is to study the ratio kjv = {T'^) /{^'^ {U))^ which decreases from 
0.572 at iV = 100 to 0.485 at iV = 500; it is far from clear whether this ratio 
is tending to zero or to a nonzero value as iV — )• oo. Clearly data at much larger 
values of N would be needed to resolve this question definitively. Unfortunately, 
such data are very time-consuming to obtain, because the deterministic algorithm 
^3 for computing T and U takes a CPU time of order iV^. Regarding the reverse 
inequality, it seems clear that pu is strictly smaller than px- 

In Table 2 we report our Monte Carlo estimates for 3-dimensional self-avoiding 
walks at iV = 100, 200, 500, 800, 1000. Here we used the pivot algorithm combined 
with deterministic algorithm ^3. For all observables, we can get a reasonable 
provided that we discard the data points at the two lowest values of N . The fits to 
(T), (T^) and var(T) yield px = 1.754 it 0.001. Again, the errors here are purely 
statistical; they do not take account of systematic errors arising from corrections to 
scaling. We will make a much more careful estimate of (T) and px in Section 4.2. 
The fits to (U), {VBarrett) and (VKarp-Luby) yield pu = 1.598 ± 0.001, 1.699 ± 0.001 
and 1.630 it 0.001, respectively. Again, the agreement is mediocre; and it is far 
from clear whether pu is strictly greater than 2px — 2 1.51 (later we will see that 
2pT — 2 sa 1.53 is a better estimate). Otherwise put, the ratio kjv = {T'^) /{^'^{U)) 
decreases from sa 0.661 at N = 100 to sa 0.528 at iV = 1000; it is far from clear 
whether this ratio is tending to zero or to a nonzero value as iV — )• oo. Again, data 
at much larger values of N would be useful. 

It is worth remarking that, in both the 2-dimensional and 3-dimensional cases, 
var(T) scales the same way as (T^) but is only about 3% (resp. 1.5%) as big. 
This means that the probability distribution of T is very narrow. Otherwise put, 
while iV-step SAWs vary radically among themselves in size and shape, the over- 
lap T((jj^^\(jj^'^'^) between tvjo of them is remarkably constant. Presumably this is 
because the operation of forming "fills in the holes" in the individual 

walks: while cu^^^ and cu^^^ are fractals, cu^^^ — cu^^^ is "semi-solid" (roughly like Swiss 
cheese). 
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4 Numerical Results 



4.1 Two Dimensions 

In Table 3 we present our data for (i?g)jv, {R^g)N^ (-^m)^? {T)n,n and the pivot- 
algorithm acceptance fraction / for SAWs in dimension d = 2 in the range 100 < 
N < 80000. Most of these SAWs were generated using the pivot algorithm, using 
either dimerization or straight rods for initialization (see Appendix B for a discussion 
of the adequacy of thermalization); run lengths were between 2 X 10^ and 8 X 
10^ pivots subsequent to thermalization. However, some data at iV < 500 were 
generated by pure dimerization (between 10^ and 2 X 10^ independent pairs of SAWs 
per run). The overlap T((jJ^^\(jJ^'^'^) was in most cases estimated using the Karp-Luby 
algorithm (Section 3.7) with 100 < -R < 200; however, some runs at iV < 500 used 
deterministic algorithm ^3 (Section 3.3). Some subtleties concerning the correct 
determination of error bars on data generated by the pivot algorithm are discussed 
in Appendix C. 

Early versions of our program had a bug, in which SAWs were pivoted only at 
sites A; > (i.e. never at the starting point). This is harmless for single SAWs, 
thanks to the lattice symmetries; but for pairs of SAWs it means that the relative 
orientation of the initial steps of the two SAWs is never altered by the algorithm. 
This causes a slight bias in the estimates of (T), especially for small N. However, 
we believe that this bias is completely negligible (i.e. much less than our statistical 
errors) for the values of N treated here; we have been unable to detect any sys- 
tematic differences between runs with and without this bug. Therefore, instead of 
throwing away the tainted data, we have simply indicated by an " in Table 3 those 
values of N for which some (not necessarily all) of the data suffers from the bug. 
Caveat lector. 

Table 4 shows the resulting values for the universal amplitude ratios {R^} / (Rl)} 
(R^) / (Rl) and The error bars are here determined using the triangle inequality; 
they are probably overestimates by a factor of 3. 

Log-log graphs of (-Rg), (-R^), (-Rm)? (^) / versus N are so straight that there 
is nothing to be gained by reproducing them here. We fit each of these quantities to 
the Ansatz AN^""^^^ [where power = 2iy for the three squared radii, 2A4 — 7 for (T), 
and —p for /] by performing weighted least-squares regressions of their logarithms 
against logiV, using the a priori error bars on the raw data points (Table 3) to 
determine both the weights and the error bars [126, Chapter 3]. As a precaution 
against corrections to scaling, we performed the fit with a lower cutoff N > Nmim 
and we tried all possible values of Nmin- The value (sum of squares of normalized 
deviations from the regression line) can serve as a test of goodness of fit. Let us 
define the "significance level" as the probability that would exceed the observed 
value, assuming the correctness of the power-law model and of the raw-data error 
bars; this probability can be read off the distribution with T> = n — 2 degrees of 
freedom (DF), where n is the number of data points used in the fit. An abnormally 
large value of (say, a level less than 5%) may indicate either that the pure power- 
law Ansatz is incorrect (e.g. due to corrections to scaling) or else that the claimed 
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error bars on the raw data are too small; further investigation would be necessary 
to determine which of these is the true cause. ^® An abnormally small value of 
(say, a level greater than 95%) probably indicates that the claimed error bars on 
the raw data are too large. 

The exponents estimated from fits to (-Rg), (-R^), (-Rm) (^) plotted 

as a function of Nmin in Figure 1. The exponent estimated from the fit to / is 
plotted in Figure 2. Looking at the values in these fits, we find statistically 
significant corrections to scaling in (-Rg) and (-R^) only for Nmin ^ 200, in (T) only 
for Nrain ^ 400, and in / only for Nmin ^ 6000. Using in each case the next larger 
value of Njnim we find: 

{Rl): V = 0.74967 ± 0.00011 
A = 0.77582 ±0.00138 

iV„™ = 300 

X^ = 32.01 (34 DF, level 

{Rl): V = 0.74963 ± 0.00008 
A = 0.10890 ±0.00014 

N.rara = 300 

X^ = 29.01 (34 DF, level 




(Ri): V = 0.74972 ± 0.00021 
A = 0.34073 ±0.00123 
iV„,„ = 1000 

X^ = 24.72 (26 DF, level = 54%) 



(T): (2A4 - 7)/2 = 0.74955 ± 0.00010 
A = 1.36283 ±0.00213 

iV^™ = 500 

X^ = 26.76 (32 DF, level = 73%) 



/: p = 0.19075 ± 0.00046 
A = 0.94095 ±0.00419 
N^,n = 7200 

= 8.29 (7 DF, level = 31%) 

(error bars are one standard deviation). The results are in excellent agreement with 
the believed exact value v = 3/4 [108,109] — although some very slight corrections 
to scaling clearly remain — and with the hyperscaling relation = 2A4 — 7. 

For (-Rg), (-Rp), (-Rm) i^)? better estimates of the (nonuniversal) amplitude 
A can be obtained by imposing the exponents 2iy = 2A4 — 7 = 3/2 and simply 

"'^^Note also that an abnormally large change in as Nmin is increased by one step — that is, a 
drop in by an amount ^1 — signals that the data point in question differs from the regression 
line by several standard deviations. This could indicate that the corrections to scaling at this value 
of N are significant, even though the overall — which is dominated by contributions from larger 
N — may look reasonable. 
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fitting observable/iV"^/^ to a constant. For all four observables, this ratio declines 
very slightly with iV; by iV 5000 it is within error bars of its apparent asymptotic 
value. Taking Nmin = 5000 we obtain: 



{Riy- 


A 


= 0.77106 ± 0.00034 








= 20.58 (24 DF, level 


= 66%) 




A 


= 0.10817 ± 0.00005 




x' 


= 20.75 (24 DF, level 


= 65%) 




A 


= 0.33894 ± 0.00015 






x' 


= 21.56 (22 DF, level 


= 49%) 



(T): A = 1.35158 ± 0.00047 

= 22.26 (24 DF, level = 56%) 

Another approach is to fit observable/iV"^/^ to a ± 6(iV/1000)-^ with a, h and A aU 
variable. For (i?^) (for which we have data only at iV > 1000), we are unable to 
obtain a decent fit with any value of Nmin available to us; the corrections to scaling 
at iV > 1000 are too weak. For (-Rg), (-R^) and (T) we obtain reasonable fits even 
for iV„,„ = 100: 

{Rl): a = 0.77100 ± 0.00040 
h = 0.00098 ±0.00050 
A = 0.843 ±0.222 

= 35.13 (35 DF, level = 46%) 

a = 0.10815 ± 0.00005 
b = 0.00011 ±0.00004 
A = 0.943 ±0.148 

= 28.04 (35 DF, level = 79%) 

(T): a = 1.35120 ± 0.00049 
b = 0.00269 ±0.00038 
A = 0.919 ±0.057 
X^ = 29.13 (35 DF, level = 75%) 

Of course, the exponents A produced by the above fits should not be taken too se- 
riously; they could well be phenomenological "effective" exponents that summarize 
the combined effects of two or more correction-to-scaling terms (e.g. the leading 
non-analytic correction N~^^ plus the analytic correction iV~^) over some particu- 
lar range of N . These fits are nevertheless useful in providing simple "interpolation 



"'^^The purpose of writing here #/1000 instead of N is to reduce the correlation between the 
estimates of h and A. Of course, 1000 can equally well be replaced by any reasonable value which is 
roughly in the middle of the N values represented by the data points. 
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formulae" that summarize our data within error bars: 

{Rl) = O.TTlOOiV^/^ + 0.33168iV°-^" (4.1) 
{R]) = O.lOSlSiV^/^ + 0.07653iV°-^" (4.2) 
(T) = 1.35120iV3/2 + 1.53518iV°-^^^ (4.3) 

In particular, we can use these formulae to compare our data with other workers' 
data at different values of N (but only for N > 100). 

The upshot is that the corrections to scaling are quite weak in the two-dimensional 
self-avoiding walk; we can barely see them at our level of precision. A serious 
study of corrections to scaling in this model would therefore require much higher 
statistics than are available here, at large but not-too-large values of N (say, 
100 ^ iV ^ 2000). We have begun such a study, and will report the results sepa- 
rately [127]. 

The corrections to scaling on the amplitude ratios (R^) / (Rl) and \I/ 

are even weaker than on the original observables (this is because the original correc- 
tions to scaling all have the same sign). The two ratios of radii have no statistically 
significant corrections to scaling, at our level of precision, even if one considers the 
true error bars to be 1/3 of those reported in Table 4. The interpenetration ratio 
does have noticeable corrections to scaling at iV < 500, but these corrections are 
very small: see Figure 3. Fitting the data for N > Nmin to a constant, we obtain 

{Rl)/{Rl) = 0.140264 ± 0.000073 
iV™. = 100 

= 2.86 (37 DF, level ^ 100%) 

= 0.439605 ± 0.000338 
iV„,„ = 1000 

= 0.74 (27 DF, level ^ 100%) 



^* = 0.66296 ±0.00043 
iV„,„ = 1000 

= 1.45 (32 DF, level ^ 100%) 

The theorists' interpenetration ratio [defined in (2.16)] is = 0.5579 ± 0.0006. 



4.2 Three Dimensions 

In Table 5 we present our data for (i?g)jv, (-R^liv, {T)n,n and the pivot-algorithm 
acceptance fraction / for SAWs in dimension d = 3 in the range 100 < iV < 
80000. Most of these SAWs were generated using the pivot algorithm, using either 
dimerization or straight rods for initialization (see Appendix B for a discussion of the 
adequacy of thermalization); run lengths were between 2 X 10^ and 4 X 10'' pivots 
subsequent to thermalization. However, some data at iV < 600 were generated 
by pure dimerization (between 10^ and 2 X 10^ independent pairs of SAWs per 
run). The overlap T((jJ^^\(jJ^'^'^) was in most cases estimated using the Karp-Luby 
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algorithm (Section 3.7) with 100 < -R < 500; however, some runs at iV < 1000 used 
deterministic algorithm ^3 (Section 3.3). Data tainted by the bug concerning (T) 
are again indicated by " (see Section 4.1). 

Table 6 shows the results for the universal amplitude ratios and 
The error bars are determined using the triangle inequality; they are probably 
overestimates by a factor of 3. 

We began by fitting (-Rg), (-R^), {T) and / to the pure power-law Ansatz AN^"'^^'' . 
The exponents estimated from fits to (-Rg), (-R^) and (T) are plotted as a function 
of Njnin in Figure 4. The exponent estimated from the fit to / is plotted in Figure 
5. Very strong corrections to scaling are apparent for all these observables; the 
exponent estimates do not appear to stabilize until one takes Nmin ^ 10'* or more. 
Nevertheless, the estimates of v do appear to be converging to a common value 
v ~ 0.5876; in particular, hyperscaling appears to be satisfied. ^° 

Another view of these same fits is presented in Figure 6: here we plot the 
estimate of u versus N~'^^ . The idea here is that the correction-to-scaling exponent 
Ai is predicted by RG calculations [104,105,106,90,91,107] to be in the vicinity of 
0.5 (see Section 5.1.2); if this prediction is correct, then each set of estimates should 
fall roughly on a straight line, at least asymptotically as Nmin oo- The data are 
consistent with this prediction, but the large fluctuations make clear that it will be 
difficult to get accurate estimates of Ai by Monte Carlo. 

Next we tried fltting (i?^), (i?2) and (T) to the Ansatz ANP"'"^'' + BNP"'"^''-^, 
where power = 2iy for the squared radii and 3iy for (T) (we now assume hyperscaling). 
For each choice of flxed exponents v and A, we determine A and B by weighted 
least-squares and record the resulting x^; we then ask which pairs (i^. A) lead to 
an acceptable x^. (Here "acceptable" is taken to mean a signiflcance level > 32% 
for T> = n — 2 degrees of freedom, where n is the number of data points. This 
corresponds roughly to confldence limits of one standard deviation.) The results 
for (-R^) and (T) are shown in Figure 7. (The result for (-Rg) is similar to that for 
(-R^), but the band of allowed values is wider.) For Nmin = 100 there are no pairs 
(i^. A) that are satisfactory for both (-R^) and (T). This means that at such small 
N one needs more than a single correction-to-scaling term to flt both observables in 
a way compatible with hyperscaling and universality. On the other hand, for larger 
Njnin one obtains a swath of acceptable pairs (i^. A), all contained in the range 
0.58715 <iy< 0.5882 and 0.36 < A < 0.80. From this analysis it seems difflcuh to 
obtain much precision on A, but it does suggest that 

ly = 0.5877 ±0.0006 (4.4) 

(subjective 68% confldence limits). 

Now suppose we impose v = 0.5877 and flt observable/iV""'^'' [where power = 
2iy for the squared radii and 3iy for (T)] to a + 6(iV/1000)-^ with a, b and A aU 
variable. For Nmin ^ 200, the estimates of A from these three flts are not consistent 

^"^Note that the rigorous inequality dv > 2A4 — 7 [see Theorem A.l and equations (A.20)-(A.22) 
in Appendix A] means that the limiting value of the upper pair of curves must be above or equal to 
the limiting value of the lower curve. 
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(this might be guessed from Figure 7). However, for Nmin ^ 300 they are consistent, 
and yield A 0.54-0.58. 

If we now return to the Ansatz ANP"""^'' + ENP"""^''-^, and take u = 0.5877 and 
A = 0.56, we obtain reasonable fits in (Rl) for Nmin ^ 150, in (-R^) for Nmin ^ 400, 
and in (T) for Nmin ^ 600. Taking Nmin = 600 we obtain: 

{Rl): A = 1.21667 ±0.00050 
B = -0.48314 ± 0.03949 

= 13.73 (27 DF, level = 98%) 

{Rl): A = 0.19455 ±0.00007 
B = -0.11432 ± 0.00465 

= 13.63 (27 DF, level = 98%) 

(T): A = 0.94477 ±0.00028 
B = 0.61687 ±0.01413 

= 12.43 (27 DF, level = 99%) 

At the very least, these fits provide simple interpolation formulae that summarize 
our data within error bars (but for N > 400 only): 

{Rl) = 1.21667iV^-^^^^- 0.48314iV°-^^^^ (4.5) 
{R^) = 0.19455iV^-^^^^- 0.11432iV°-^^^^ (4.6) 
(T) = 0.94477iV^-^'^^^ ± 0.61687iV^-2°^^ (4.7) 

If we take these fits seriously, we obtain the results 

4\V4'] = 1.48 ±0.18 (4.8) 

= -1.64 ±0.17 (4.9) 

for the universal ratios of correction-to-scaling amplitudes. 

Another way to study the corrections to scaling is to look at an amplitude 
ratio which is known to tend to a nonzero constant as iV — )• oo, and fit it to 
a ± 6(iV/1000)~^ with a, b and A all variable. Two candidates are {R^) / {Rl) and 
(assuming hyperscaling) \E'. Unfortunately, the corrections to scaling in {Rl^) / {Rl) 
are too weak to yield much information on A (see Table 7): the amplitude b is 
within 2-3(7 of zero, and the error bars on A are enormous. However, we can obtain 
a reasonable estimate of the universal limiting value (-R^)/ {Rl) = a. Note first that 
the unusually small confirms our belief that the true error bars on (i?^)jv/ {RI)n 
are roughly 1/3 of those reported in Table 4. If so, then the error bars on a, b and 
A in Table 7 should also be reduced by a factor of 3. Making this adjustment, 
we conclude that 

(Rl) 



(Rl) 



0.1599 ±0.0001 ± 0.0001 ; (4.10) 
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here the first error bar represents systematic error due to uncontrolled corrections 
to scaling (subjective 68% confidence limits) and the second error bar represents 
the (adjusted) statistical error (classical 68% confidence limits). 

The analysis of \I/ is much more favorable for estimating A, as the corrections 
to scaling are rather strong (Figure 8). [This is because the corrections to scaling 
on (T) and (-R^) have opposite signs, as is evident from Figure 4.] A sample of 
the results is shown in Table 8. Note that the estimates of b and A are quite 
stable as Nmin is increased, and that the error bar on A is remarkably small. (The 
unusually small confirms our belief that the true error bars on \I/ are roughly 
~ 1/3 of those reported in Table 4. If so, then the error bars on 6 and A should 
also be reduced by a factor of 3.) The excellent fit can be seen graphically in 
Figure 9; \I/ is amazingly close to a linear function of iV~°'^^. This confirms our 
expectation that the correction-to-scaling exponent Ai is approximately 0.5, and 
suggests that it may be slightly higher, around 0.56. A fair estimate would be 

Ai = 0.56 ±0.03 (4.11) 

(subjective 68% confidence limits). On the other hand, the A produced by the 
above fit could well be a phenomenological "effective" exponent that summarizes 
the combined effects of two or more correction-to-scaling terms (e.g. the leading 
correction N~^^ plus the analytic correction iV~^) over some particular range of N . 
The only way to sort this out would be to use a larger Nmin together with improved 
statistics. In any case, we can estimate the universal limiting value 

^* = 0.2471 ± 0.0001 ± 0.0002 . (4.12) 

(We have again made the factor-of-3 adjustment in the statistical error bar.) The 
theorists' interpenetration ratio [defined in (2.16)] is = 0.2322 ± 0.0001 ± 0.0003. 
Since \E'* > 0, hyperscaling is satisfied. 

5 Discussion 

5.1 Comparison with Previous Numerical Studies 

5.1.1 Two Dimensions 

We first compared our raw data with those of other studies [128,10,129,130,127], 
making direct comparisons where the values of N match and using the interpolation 
formulae (4.1)-(4.2) at other values of iV.^^ We find excellent agreement, with one 
perplexing exception: when compared with the extremely precise data of [127], our 
values for (Rl) and (-R^) are about 3.5cr low at iV = 100, and slightly low (< 2a) 
at 200 < iV < 500. Interestingly, these are precisely the runs that were performed 

^^Warning: We believe that the theoretical premises of [130] are erroneous, in that the authors of 
[130] fail to distinguish correctly which quantities are universal and which are non-universal. (See 
Section 1.2 above.) Nevertheless, the Monte Carlo data in [130] are useful. 



35 



using pure dimerization, while the runs at iV > 500 used the pivot algorithm. It is 
therefore conceivable that our dimerization program has a subtle bug which causes 
a very small error (perhaps one that decreases with N). However, we have been 
unable to find any such bug, and the small discrepancy could well be a statistical 
fiuke. 

Our estimates of the universal ratios and agree perfectly 

with the best published estimates [129], and have roughly the same precision. In 
particular, they confirm the beautiful conformal-invariance prediction of Cardy and 
Saleur [131] (as corrected by [129]) 

91 im {RD + 2 - ° • ^'-'^ 

Much more precise data will be available soon [127]. 

To our knowledge there are no previous series-extrapolation or Monte Carlo 
estimates of \E'* for two-dimensional polymers. But we can compare our value 
with field-theoretic estimates: 

Monte Carlo, simple cubic lattice, 1000 <N< 80000 (this work): 0.66296 ± 0.00043 

Edwards model through order e^, naive sum (Appendix D below): 0.8243 
Edwards model through order e^, exploiting d = 1 value (Appendix D below): 0.6647 

This last estimate is amazingly close to the correct value; it would be interesting to 
know whether this close agreement is an accident. 



5.1.2 Three Dimensions 

We first compared our raw data with those of other studies [133,10,130,134,135,116,136,137], 
making direct comparisons where the values of N match and using the interpola- 
tion formulae (4.5)-(4.6) at other values of N We find excellent agreement. We 
also compared our raw data with Nickel's interpolation/extrapolation formulae [13, 
equations 6 and 7] , which are based on exact enumeration of short chains combined 
with high-precision Monte Carlo data at iV < 2560. The agreement for (-Rg) and 
(i?^) is excellent: most of the points differ by less than Icr, about 20% differ by 
between Icr and 2cr, and none differ by more than 2a. However, for (T) there are 
some modest discrepancies: for N < 600 and N = 800, our values are all between 
2(7 and 5.2cr lower than Nickel's; while for N > 15000, about half of our values are 
between la and 3a higher than Nickel's. We do not know whether either of these 
discrepancies is real, or what might be its cause. We suspect that it arises simply 
from statistical errors in Nickel's raw data for (T) (which might be several times as 
large as our errors); these would induce statistical errors in the coefficients of his 

^^Table 9 (see Section 5.2 below) gives the values of for N < 7, from exact enumeration [132]. 
But this series is much too short to be usefully extrapolated. 

^^Warning: Reference [116] uses an unconventional definition of (i?^); it is (1 -|- l/N)^ times ours. 
Note also that the data of references [134,135,137] lack error bars. 
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extrapolation formula, which would in turn induce correlated statistical errors at 
nearby values of N . 

Next we compare our estimates of the critical exponents v and Ai with previous 
work: 

Monte Carlo, simple cubic lattice, 100 < iV < 80000 (this work): 
// = 0.5877 ± 0.0006, Ai = 0.56 ± 0.03 

Monte Carlo, simple cubic and BCC lattices, 120 < iV < 2400 [133]: 
// = 0.592 ±0.002 

Monte Carlo, simple cubic lattice, 200 < iV < 3000 [10]: 
// = 0.592 ±0.002 

Monte Carlo, simple cubic lattice, 200 < iV < 7168 [116]: 
// = 0.5909 ± 0.0003 

Series extrapolation, various lattices [138]: 
// = 0.592 ±0.003 

RG, ra = field theory [104,105,106]: 

// = 0.5880 ± 0.0015, Ai = 0.47 ±0.025 

RG, ra = field theory [lOl]'^'^: 

V = 0.5872 ± 0.07(5* - 1-39) ± 0.0004 

Edwards model through order [90]^^: 
// ^ 0.588, Ai ^ 0.473 

Edwards model through order [91]^^: 
z/ ^ 0.5886 Ai ^ 0.465 

(all Monte Carlo error bars are one standard deviation). It is clear in retrospect 
that the earlier Monte Carlo estimates [133,10,116], based on shorter walks than 
those used here, were biased upwards due to corrections to scaling. (This effect is 
seen very clearly in Figure 4. Moreover, if we truncate our own data to lie in the 
same range of N as the previous studies, the resulting estimates of v are almost 
identical to the quoted ones.) We have now done what we believe is a careful analysis 
of the corrections to scaling, and we have obtained reasonably good control over 
them. We therefore think that the current estimate is correct within its claimed 
error bar. (But we could be wrong!) We do not know why the series-extrapolation 
estimates are also high, but it could arise from the same effect; perhaps the walks 
probed in these analyses (up to 20-30 steps) are simply not long enough to permit 
an adequate analysis of corrections to scaling, even using the most sophisticated 
differential- approximants methods. 

^■*The best estimate of [107] for g* is 1.39. Comparison of references [104] and [107] suggests that 
the uncertainty in g* may be of order 0.03-0.04. This would add an extra uncertainty of order 
0.002-0.003 to v. 

^®The term of order has recently been obtained [107], but the extended series has not yet (to 
our knowledge) been analyzed. 
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The amazing fact is that the Monte Carlo estimates of v have stabihzed at almost 
exactly the value predicted by the field-theoretic calculations in their various equiv- 
alent forms (n = field theory [104,105,106,107] or Edwards model [89,90,91,107]). 
The very high accuracy of the field-theoretic calculations is rather surprising to us, 
since this method is susceptible to serious (and quite possibly undetectable) sys- 
tematic errors arising from a confluent singularity at the RG flxed point [33,139]. 
The weakness of this effect may be related to the apparent fact [140,139] that the 
confluent exponents 1/Ai and A2/A1 are both very close to an integer (namely, 
they are 2). 

On the other hand, for the confluent exponent Ai, the agreement between our 
Monte Carlo estimates and the fleld-theoretic predictions is not so good. Our Monte 
Carlo data for ^ can be flt amazingly well, all the way down to iV = 100, by the 
Ansatz ^ = a+h/N^ with A = 0.561±0.004 (statistical errors only). This exponent 
differs considerably from the fleld-theoretic predictions, which all lie at Ai 0.47. 
We do not know which of these estimates, if either, is the correct one — indeed, 
there are good reasons to be distrustful of both! On the one hand, our Monte Carlo 
estimate could well be a phenomenological "effective" exponent that summarizes 
the combined effects of tvjo or more correction-to-scaling terms (e.g. the leading 
correction N~^^ plus the analytic correction iV~^) over some particular range of 
N . (To test this, we would need to go to larger Nmin and obtain signiflcantly 
improved statistics.) As things stand, our data do not rule out the fleld-theoretic 
prediction Ai 0.47, provided that one includes suitable subleading correction- 
to-scaling terms. On the other hand, the fleld-theoretic prediction for Ai should 
also be taken with a grain of salt: it arises from the slope of the /3-function at the 
flxed point g = </*, and as Nickel [33,139] showed long ago, this function has, in 
addition to the desired term linear in g — g* , also nonanalytic terms proportional to 
{g* — gY^^^ , (g* — g)^^^^^ and the like (see Section 5.2 below for further discussion). 
But the numerical methods currently employed [104,105,106,107,139,89,90,91] to 
extrapolate the perturbation series from g = to g = g* assume, contrary to fact, 
that l3(g) is regular at g = g* . The presence of terms (g* — gY with 1 < p < 2 could 
well lead to systematic errors in estimates of the slope fi'{g*) [33]. It is therefore 
worth considering the possibility that the true Ai is indeed closer to sa 0.56 than 
to sa 0.47, and that the error bars on the fleld-theoretic estimates may be over- 
optimistic. 

Let us also make a brief comparison with the experimental results on polymers 
in a good solvent; a more detailed analysis will appear elsewhere [141]. The most 
systematic data on the static scaling behavior of high-molecular- weight polymers in 
a good solvent were obtained in the 1970's by three Japanese groups [142,143,144], 
using polystyrene in benzene. These data were reanalyzed, shortly after the appear- 
ance of the RG predictions for v [145], by Cotton [146]. After making corrections 
for polydispersity. Cotton obtains the value v = 0.586 it 0.004 — in good agreement 
with the RG prediction, and now in good agreement with the Monte Carlo data as 
well. Unfortunately, we believe that this experimental value is unreliable, for the 
following reasons [141]: 
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1) The samples of Yamamoto et al. [142] have an unknown polydispersity, but 
this polydispersity is certainly not zero (as Cotton's analysis implicitly as- 
sumed). 

2) The measurements of Fukuda eA al. [143] were apparently afflicted by a serious 
systematic error (of magnitude ranging from 5% to 18%) arising from the way 
that the solution concentration was measured: see [147, footnote 19].^^ 

3) The data from the three different laboratories cover almost-nonoverlapping 
ranges of molecular weight. As a result, the combined analysis of the three 
sets of data is highly susceptible to the effects of small systematic discrepancies 
in absolute calibration between the three laboratories, as well as to the fact 
that the three laboratories used slightly different temperatures (30° C for 
[142,143] and 25° C for [144]). 

4) The raw data lack error bars. As a result, it is impossible to distinguish cor- 
rections to scaling (or systematic experimental errors) from statistical errors. 
This distinction is crucial to extracting a reliable value for as our analysis 
of Monte Carlo data has demonstrated. 

A new generation of experiments, using modern ultra-sensitive light-scattering in- 
strumentation [149,150,151] and an optimal statistical analysis [152], could prove 
to be very exciting. 

Next let us compare our estimates of the universal ratio (-R^)/ (-Rg) with previous 
work: 

Monte Carlo, simple cubic lattice, 100 < iV < 80000 (this work): 0.1599 ± 0.0001 ± 0.0001 

Monte Carlo, simple cubic lattice, 120 < iV < 2400 [133]: 0.1597 ± 0.0003 

Monte Carlo, body-centered-cubic lattice, 120 < iV < 2400 [133]: 0.1594 ± 0.00015 

Monte Carlo, simple cubic lattice, 200 < iV < 3000 [10]: 0.1603 ± 0.0004 

Monte Carlo, simple cubic lattice, 200 < iV < 7168 [116]: 0.1596 ± 0.0002 

Monte Carlo, simple cubic lattice, 26 < iV < 3328 [136]: 0.16003 ± 0.00003 

Edwards model through order z'^ [136]: 0.16012 ± 0.00030 



(all Monte Carlo error bars are one standard deviation). The agreement is excellent.^ 
It is worth noting that is non-monotonic as a function of N: for small 

N (^ 15), exact enumeration [154,137] shows that {R'^) / (Rl) is considerably above 
its limiting value 0.1599; but for larger N (e.g. 100 ^ iV < 1000), our data in 
Table 6 show unequivocally that (R^) / (Rl) is below its limiting value; and our ffl 



^^Unfortunately, it will not suffice simply to replace the measurements of Fukuda ei al. [143] by 
those of Utiyama ei al. [147] on a subset of their samples, because these latter measurements may 
suffer from a systematic error of their own, arising from the extrapolation to zero scattering angle: 
see [148, footnote 30]. 

^^There also exists a calculation in the Edwards model to second order in e = 4 — c? [153], but 
unfortunately this expansion is too ill-behaved to be extrapolated reliably to e = 1 (see [10, p. 183]). 
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(Table 7) suggests that as iV — )• oo the approach to the hmiting value is also from 
below. It would be interesting to know whether series-extrapolation techniques can 
"sense" this nonmonotonicity (which lies beyond the present enumerations) and 
predict approximately the correct limiting value. 

Next let us compare our estimates of the limiting interpenetration ratio with 
previous work: 

Monte Carlo, simple cubic lattice, 100 < iV < 80000 (this work): 0.2471 ± 0.0001 ± 0.0002 

Monte Carlo, simple cubic lattice, 13 < iV < 2560 [13,155]: 0.2465 ± 0.0012 

Edwards model through order z'^, Pade [1/1] (Appendix D below): 0.2092 
Edwards model through order e^, naive sum (Appendix D below): 0.2686 
Edwards model through order e^, exploiting d = 1 (Appendix D below): 0.2486 



(all Monte Carlo error bars are one standard deviation).^® Our estimate is thus 
consistent with the very recent Monte Carlo estimate of Nickel et al. [13,155], but 
is about 4 times as precise. The Edwards-model (renormalization-group) estimates 
are all of the right order of magnitude, but most of them are not terribly close to the 
correct answer. (This is hardly surprising, in view of the very short perturbation 
series on which these estimates are based.) The e-expansion result augmented by 
exact information at c/ = 1 is, however, amazingly close to the correct value, both 
in (/ = 3 and d = 2. It would be useful to obtain a better understanding of whether 
this is a coincidence or not — perhaps by calculating the Oic") term in \E'*. 

Finally, it is worth noting that the experimental values for \E'* also lie in the range 
0.22-0.25 [11, Section 10. F] [96, Sections 2.2.5 and 2.2.6] — or more optimistically, 
they average to 0.245 it 0.005 [146]. However, the experimental measurements of 
\E' are subject to all the problems noted earlier, as well as the danger of addi- 
tional serious systematic errors arising from curvature in the extrapolation to zero 
concentration. This is a particularly severe problem for older studies which used 
higher concentrations in order to compensate for the less-sensitive light-scattering 
instrumentation then available. Again, new experiments would be highly desirable. 

As explained in Section 1.2, also the ratios of correction-to-scaling amplitudes 
^rI/^rI ^A^ /^rI universal. So ideally we would like to measure these am- 
plitude ratios, and compare them with the Edwards-model prediction h^R^/h^^^ = 
1.249 lb 0.035 [136]. Unfortunately, the correction-to-scaling amplitudes are ex- 
tremely sensitive to the choice of exponents v and Ai. Therefore, the only sensible 
comparison that can be made is to impose the same values of v and Ai as are 
used in [136] — namely, v = 0.588 and Ai = 0.47 — and then perform a weighted 
least-squares fit to estimate the amplitudes. Unfortunately, the resulting estimates 
of the correction-to-scaling amplitudes are highly unstable as a function of Nmin- If 
we take iV„,„ = 3000, we obtain b^^] = -0.10 ± 0.05 and = -0.17 ± 0.05, hence 

^R^ /^rI — '^•^ i l-S — ^ very useful estimate! 



^®Note also that Table 9 (see Section 5.2 below) gives the values of for N < 7, from exact 
enumeration [132]. But this series is much too short to be usefully extrapolated. 
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Finally, we can compare the estimates of the pivot-algorithm acceptance-fraction 
exponent p: 

Monte Carlo, simple cubic lattice, 100 < N < 80000 (this work): p ~ 0.105 or perhaps lower 

Monte Carlo, simple cubic lattice, 200 < iV < 3000 [10]: p = 0.1069 ± 0.0009 

Monte Carlo, tetrahedral lattice, 49 < iV < 9999 [115]: p ^ 0.1146 

Monte Carlo, simple cubic lattice, 100 < iV < 1600 [135]: p ^ 0.103 

There is a slight discrepancy between our results and the very careful work of Zifferer 
[115]; we do not understand its origin. It is conceivable (though in our opinion un- 
likely) that the critical exponent for the pivot-algorithm acceptance fraction might 
vary from one three-dimensional lattice to another. More likely, the magnitude (and 
even the sign) of the corrections to scaling may vary radically between lattices. 

5.2 The Sign of Approach to ^* 

For several decades, most work on the behavior of long-chain polymer molecules 
in dilute solution [156,157,158,11,12,96] has been based on the "two-parameter the- 
ory" in one or another of its variants: traditional (Flory-type)^^, pseudo-traditional 
(modified Flory-type)"^° or modern (continuous-chain-type)"^^. All two-parameter 
theories predict that in the limit of zero concentration, the mean-square end-to-end 
distance (-Rg), the mean-square radius of gyration (-R^) and the interpenetration ra- 
tio \I/ depend on the degree of polymerization N (or equivalently on the molecular 
weight M = N M monomer ) accordiug to 

{Rl) = ANFRjbN) (5.2a) 
{RD = ANFn,{bN) (5.2b) 
^ = F^{bN) (5.2c) 

where Fr^^ Fr^^ Fq, are claimed to be universal functions (which each specific two- 
parameter theory should predict), and A and b are non-universal scale factors de- 
pending on the polymer, solvent and temperature but independent of N. [The 
conventional notation is a|j = -Fr^, a| = 6-FRg, h = agFi^, / z and z = (6iV)^~''/^ in 
spatial dimension d.] Moreover, virtually all the theories — and in particular the 
modern continuous-chain-based theories — predict that Fq, is a monotone increas- 
ing and concave function of its argument 6iV, which approaches a limiting value 
^* ?^ 0.2 - 0.3 (for d = 3) as bN ^ oo. 

^^See Yamakawa [157], Sections 11 and 16 (pp. 69-73 and 110-118) and parts of Sections 15, 20b 
and 21b (pp. 94-110, 153-164 and 167-169). See also DesCloizeaux and Jannink [12], Section 8.1 
(pp. 289-313). 

^"^See Yamakawa [157], most of Section 15 (pp. 94-110) and parts of Sections 20b and 21b (pp. 153- 
164 and 167-169). See also Domb and Barrett [159]. 

^^These theories take as their starting point the Edwards model of a weakly self-avoiding contin- 
uous chain [83,84,85,86,87,11,12]. (The Edwards model is also equivalent to the continuum ip^ field 
theory with n = components.) See DesCloizeaux and Jannink [12] for a detailed treatment of the 
Edwards model. 
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But our Monte Carlo data for the self-avoiding walk (Figures 3 and 8) show 
precisely the opposite behavior: \I/ is a decreasing and convex function of iV, which 
approaches a hmiting value 0.2471 {d = 3) or 0.6630 {d = 2) as N ^ oo. 

The same behavior was found by Nickel [13]. The decrease of \I/ with N is strong 
in (/ = 3, and weak (but noticeable) in d = 2. 

In retrospect, this behavior is heuristically almost obvious: Short self-avoiding 
walks behave roughly like hard spheres, i.e. \I/ is on the same order of magnitude 
as the hard-sphere value (2.20) [see Table 9]. On the other hand, long self-avoiding 
walks are fractal objects, i.e. "thinner" than hard spheres, so one expects <^ 
hard- sphere- This is manifestly true in dimension d > 4 (where = 0) and in 
dimension d = 4 — e [where = e/8 -|- 0(e^)]; it is natural to expect (and we now 
confirm) that it is true also in c/ = 3, and to a lesser extent in d = 2. Of course, 
in (/ = 1, SAWs are hard spheres (i.e. hard rods), so = hard- sphere ■^'^ (The 
monotonic decrease of /"^ hard-sphere as a function of d is shown in the last line 
of Table 9.) If one now conjectures the simplest behavior, namely that \E'jv is a 
monotonic function of iV, it follows that (in dimension d > 2) \E'jv must approach 
its limiting value \E'* from above. 

There is also experimental evidence [93,94,95,96] that for real polymers in a 
sufficiently good solvent, the approach to \E'* is from above, contrary to the two- 
parameter theory. This behavior was considered to be a perplexing "anomalous 
effect", and various purported explanations were advanced [160,161,94]. 

The correct explanation, in our opinion, was given three years ago by Nickel 
[13] (see also [14,15]): theories of two-parameter type are simply wrong. Indeed, 
they are wrong not merely because they make incorrect predictions, but for a more 
fundamental reason: they purport to make universal predictions for quantities that 
are not in fact universal. Two-parameter theories"^"^ predict, among other things, 
that \E' is a universal function of the expansion factor a| = {R^} / {Rl)Tg] in particu- 
lar, \E' is claimed to depend on molecular weight and temperature only through the 
particular combination a|(M, T). This prediction is quite simply incorrect, both 
for model systems and for real polymers. Indeed, even the sign of the deviation 
from the limiting value \E'* is not universal. 

All this has a very simple renormalization-group explanation [13,14,15], so it is 
surprising that it was not noticed earlier. As mentioned already in Section 1.2, stan- 
dard RG arguments predict, for any real or model polymer chain, the asymptotic 
behavior 

(Rl) = A«.iV^ll + + • • •) (5.3a) 

(i?^) = An,N'^{l + b^^lN-^^+...) (5.3b) 

^ = ^* (1 + + . . .) (5.3c) 

^^In d = 1, is essentially constant with N; it differs from ^* only by corrections of order l/N 
arising from the discreteness of the chain. The exact formula is = C^/ V^)(J^ ~^ \)l \/^(J^ + 2)- 

^^More precisely, two-parameter theories in which the scale factor A is independent of temperature. 
In [15] these are termed "strong two-parameter theories". 
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as iV — )• oo at fixed temperature T > Tg. The critical exponents v and Ai are 
universal. The amplitudes Ar^, ^Rgi ^rI^ ^rI^ nonuniversal; in fact, even the 

signs of the correction-to-scaling amplitudes 6^], b^^^ are nonuniversal. However, 
the RG theory also predicts that the dimensionless amplitude ratios Ar^/Ar^^ 

^rI/^rI ^^^^ /^rI universal [13,88]. 

Recently, however, several papers have appeared [162,163,164,165] which at- 
tempt to explain the observed approach to from above in terms of either an al- 
leged "second branch" of the renormalized field theory with renormalized coupling 
constant g > the fixed-point value g* [166,164,165] or an "extended two-parameter 
theory" containing an extra parameter [162,163]. We confess that we have been 
unable to understand the conceptual basis of these papers. Here is our attempt 
[14,15] to clarify what is going on: 

We believe that critical phenomena are best understood in a Wilson-type (or 
"Wilson-de Gennes-type") renormalization-group framework [167,158]. By this we 
mean an RG map TZ acting on the (infinite-dimensional) space of Hamiltonians 
for a field theory or polymer model with some (fixed) ultraviolet cutoff A < oo 
(e.g. living on a lattice); the map TZ acts by integrating out the high-momentum 
(short-wavelength) degrees of freedom. "^"^ We wish to contrast this with a "field- 
theoretic" RG [69,11,12], which acts in a (finite-dimensional) space of continuum 
(renormalized) field theories or polymer models (i.e. models with the ultraviolet 
cutoff already taken to infinity) by spatial dilation. A few paragraphs from now we 
will explain the field-theoretic RG in terms of the Wilson RG (and we will see that 
the reverse is not possible). But for the moment let us simply proceed with the 
Wilson approach. 

Figure 10 (taken from [14]) shows part of the (Wilson) renormalization-group 
flow for a cutoff fleld theory or polymer model in dimension 3 < c/ < 4. Here 
H* is the Gaussian flxed point (which for 3 < c/ < 4 is also the flxed point con- 
trolling the theta-solvent behavior), while Hq^ is the nontrivial (= good-solvent 
= Wilson-Fisher) flxed point. (Please ignore for now the curves Ais^ Mu and C.) 
More precisely. Figure 10 shows the flow on the critical surface (corresponding to 
correlation length = oo in the cutoff fleld theory, or chain length iV = oo in the 
polymer model). Non-critical models (i.e. < oo or iV < oo) lie above the plane 
of the page. Both H* and Hq^ have unstable (= relevant) directions coming out 
of the page, and these trajectories lead to the inflnite-temperature flxed point H'^ 
(which has = or iV = 0). 

We must now distinguish two very different limiting situations in polymer the- 
ory: 

(I) iV — )• oo at flxed temperature T, where either (a) T > Tg^ (b) T = Tg, or (c) 
T<Te. 



^■*For example, for a continuous-space field theory with ultraviolet cutoff A, the map TZ could be 
integration over field components with momenta in the shell A/2 < |p| < A [167]; for a lattice field 
theory, TZ could be passage to block spins [168]; for a discrete polymer chain, TZ could be decimation 
of odd-numbered sites along the chain [158]. 
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(II) N — > oo, T —i- Tg with x = N'^{T — Tg) fixed, where ^ is a suitable crossover 
exponent. 

Roughly speaking, case la corresponds to the good-solvent regime, while case II 
corresponds to the crossover regime near the theta point. "^^ 

Case la applies to any one-parameter family of Hamiltonians, parametrized by 
the chain length iV, that transversally intersects the critical surface (i.e. the plane 
of the page) at some point within the domain of attraction of Hq^ (i.e. anywhere to 
the right of Ais)- For example, the family could intersect the critical surface at 
P , Q ov R (among many other places). Then, the critical exponents and universal 
amplitude ratios [cf. (5.3) ff.] associated with the limit iV — )• oo in the family 
are completely determined by the RG flow in an inflnitesimal neighborhood of 
Hq^] they are therefore the same for all families intersecting the critical surface 
within the domain of attraction of Hq^ ("universality"). On the other hand, the 
nonuniversal amplitudes arise from the entire history of the RG flow that takes 
to Hq^] they are therefore different for different families JF. For example, a family 
crossing the critical surface at P or i? would approach \E'* from below, while a family 
crossing at Q would approach \E'* from above. The nonuniversal amplitudes cannot 
be predicted except through detailed knowledge of the microscopic physics (i.e. of 
the family JF). In particular, they cannot be predicted by any "coarse-grained" 
theory (such as a renormalized fleld theory), or indeed by any simple mathematical 
model. 

As will be explained below, the manifold (which extends also out of the 
plane of the page) corresponds to the continuum Edwards models, with the plane 
of the page corresponding io z = oo. It is important to note that has no special 
status with respect to the good-solvent flxed point Hq^] it is merely one of many 
trajectories whose intersection with the critical surface is attracted to Hq^. The 
universal properties of polymer chains in a good solvent (case la) can indeed be 
extracted from the limit 2; — )• oo of the continuum Edwards model; but they can be 
extracted equally well from the limit iV — )• 00 of any family JF whose intersection 
with the critical surface is attracted to Hq^. One family may be more convenient 
for computation than another, but all have the same conceptual status. 

Let us now explain the relation between the Wilson RG and continuum fleld 
theories. Let H* be any critical flxed point (for the moment it doesn't matter 
whether H* is Gaussian), and let M.^ (resp. M.u) be the stable (resp. unstable) 
manifold of H*.^^ Continuum limits are obtained by taking a sequence of initial 
Hamiltonians Hn approaching the stable manifold, and rescaling lengths by suitable 
factors. This rescaling is equivalent to applying the map TZ a suitable number 



^^However, it is crucial to understand that cases I and II refer to families of Umttmg paths in the 
(T, #)-plane, not to "regions" or "domains" of the finite (T, #)-plane. Failure to appreciate this 
distinction — which may at first seem rather pedantic — can lead to apparent paradoxes regarding 
the sign of the correction-to-scaling terms. 

^®For simplicity let us assume that there are no marginal operators. The presence of marginal 
operators (such as the ip^' operator at the Gaussian fixed point in c? = 3) does not affect the 
fundamental conclusions, but merely induces multiplicative logarithmic corrections. 
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of times. The low-energy effective Hamiltonians H!lf = TZ^Hn then tend to the 
unstable manifold (see Figure 10). Continuum field theories F are thus in one-to- 
one correspondence with Hamiltonians H on the unstable manifold: the correlation 
functions of F at momenta |p| < A are equal to the correlation functions of H with 
cutoff A. This point of view has been emphasized by Wilson and Kogut [167] and 
others [169,168]. 

In particular, the Gaussian fixed point H* has (in dimension 3 < c/ < 4) a 
two-dimensional unstable manifold Aiu- the two unstable directions correspond to 
(fi'^ and (fi"^ interactions (or in polymer language, a chain-length fugacity and a two- 
body self-interaction)."^'' So Aiu corresponds to the manifold of superrenormalizable 
(fi"^ continuum field theories — or in polymer language, to the continuum Edwards 
model. 

We can now understand the connection between the Wilson and field-theoretic 
renormalization groups."^® The homogeneous field-theoretic RG equations [171,172,173] 
describe how a family of continuum field theories is mapped into itself under spatial 
dilation. On the other hand, the continuum field theories are in one-to-one corre- 
spondence with Hamiltonians on the unstable manifold; and this correspondence 
takes spatial dilation into the RG map TZ. Thus, the field-theoretic RG is nothing 
other than the Wilson RG restricted to the unstable manifold and then rewritten in 
terms of "renormalized" parameters. 

Having understood this, we can now evaluate the attempts [166,164,165] to 
explain the "wrong" sign of approach to (or the analogue in liquid-gas critical 
points) in terms of an alleged "second branch" of the renormalized field theory, 
located at g > g*.^^ The trouble is that, as far as we know, no such branch exists. 
As explained above, continuum field theories correspond to the unstable manifolds of 
critical fixed points. Thus, the alleged "second branch" could exist only if there were 
an as-yet-unknown critical fixed point H*^^ lying somewhere to the right of Hq^^ 
part of whose unstable manifold is attracted to Hq^. There is (to our knowledge) 
no evidence whatsoever for the existence of such a fixed point. 

^^We ignore here the marginal operator (Vi,c)^, which corresponds to a physically trivial rescaling 
of field strengths. 

^®This connection was worked out by Kupiainen and Sokal in 1983, and was quite possibly known 
to others as well. However, to our knowledge it first appeared publicly in [14]. Very similar ideas 
appeared earlier in work of Hughes and Liu [170]. 

^^Of course, no such explanation is needed, because we have already given a complete and straight- 
forward explanation in terms of the Wilson RG. But there would be no harm in giving an alternate 
explanation of the same phenomenon — provided that this explanation is correct! 

■*'^The nonexistence of such a branch of continuum field theories seems to be recognized by Kriiger 
and Schafer, who note [164, abstract and p. 764] that "the strong coupling branch implicitly relies 
on the existence of a finite segment size" (~ ultraviolet cutoff). But on the other hand they insist 
that renormalized field theory can be used for g > g* [165, p. 3]. We do not understand how 
these two statements can be reconciled. Perhaps Kriiger and Schafer want to study the cutoff field 
theories to the right of iJ^^. But there is no distinguished one-parameter family of such theories (i.e. 
the putative extension of A4u)', there is simply the infinite-dimensional space of all cutoff theories 
(with complicated Hamiltonians as well as simple ones) lying to the right of H^g, all of which have 
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The conventional field-theoretic RG approach [69,166,164,165] makes another 
error in assuming that the /3-function describing the field-theoretic RG fiow is reg- 
ular at g = g* . In fact, as pointed out long ago by Nickel [33,139], the /3-function 
should be expected to contain nonanalytic terms like (g* — gY^^^^ (g* — g)^^^^^ and 
many others. This can easily be understood in our Wilson-type RG framework: The 
field-theoretic /3-function describes the Wilson RG fiow restricted to the unstable 
manifold Aiu (and rewritten in terms of the "renormalized" parameter g). Now Aiu 
has no special status at Hq^: like every other RG trajectory'*^, it approaches Hq^ 
tangent to the leading irrelevant direction; but barring a miraculous coincidence, 
Aiu also has nonzero components (with respect to the nonlinear scaling fields at 
Hq^) in all of the subleading irrelevant directions. This induces nonanalytic terms 
(</* ~ fl')^^^^S (</* ~ fl')^^^^S • • • in the /3-function. In addition, the analytic correc- 
tions to scaling at Hq^ induce nonanalytic terms like (g* — gY^^^^ (g* — gY^^^^ ■ ■ ■ 
in the /3-function [33]. 

In summary, the field-theoretic /3-function is nonanalytic as g '\ g*^ and is (as 
far as we know) not defined at all for g > g*. These facts cannot be understood by 
manipulation of formal expressions within renormalized field theory, but they can 
be understood when the renormalized field theory is placed within the Wilson-type 
RG context. 

Finally, this framework allows us to understand [14,15] the special role played 
by the continuum Edwards model. This model has no special status at Hq^^ but it 
does have special status at H*: it is the unstable manifold Aiu- For this reason, the 
continuum Edwards model describes the universal crossover scaling behavior in an 
infinitesimal region just above the theta temperature ^ namely the limit N — )• oo, T — > 
Te with X = N^{T - Te) fixed > (case II above), where = 2 - c//2 for 3 < c/ < 4 
and ^ = I X log""*^** for d = 3."*^ That is, the continuum Edwards model controls 
the behavior of any two-parameter family JF' of Hamiltonians that transversally 
intersects the critical surface in some curve C which in turn transversally intersects 
Ais (see Figure 10). Since the theta point is beyond the scope of the present paper, 
we refer the reader to [14,15] for details. 

equivalent conceptual status. 

The same objection applies to the "extended two-parameter theory" of Chen and Noolandi 
[162,163]. Perhaps their extra parameter is intended to correspond to the coefficient of the leading 
irrelevant coupling at H^g] but in that case it is merely an inordinately complicated restatement of 
(5.3) ff., and moreover it neglects the second and higher irrelevant couplings. Of course, one could 
consider the manifold corresponding to the vanishing of the second and higher irrelevant couplings 
(with respect to the nonlinear scaling fields at H^g); but this manifold, on the left side of H^g, 
does not coincide with A4u- 

Finally, we believe that the same objection applies to the second half of Nickel's paper [13, after 
p. 1359], in connection with his "recursion model". He seems to realize this, as he describes the 
upper branch in his simplified recursion model as "(/aasmniversal" . 

Except for a measure-zero set of exceptional trajectories. 

■*^By the latter expression we mean that in c? = 3 the correct scaling variable is x = 
#i/2(log#)-4/ii(T - Te) [174,175,176]. 
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5.3 Prospects for Future Work 



Let us close this paper by mentioning briefly some interesting areas for future 
work. 

Higher virial coefficients can be studied by the same methods that we have used 
here to study the second virial coefficient. For example, the third virial coefficient 
between molecules of types z, j, k is [102] 

i?P = I E E W,{s)W,{s')W,{s") X 
s e s, x' ei'^ 
s' e Sj x" e TL'^ 
s" e 

[l _ g-v,,((o,.),(^',.'))] _ g-v,,((o,.),(^",."))j _ g-v,,((^',.'),(^",«"))] 

(5.4) 

[compare (2.13)], where the interaction energy V is given as usual by (2.14). The 

ratio B^^'^'^^ / B2^'^^ is expected to be universal for polymers in a good solvent, 
and it would be interesting to know its value. (At present, only a crude field- 
theoretic estimate is available: B3/B2 ~ 0.517, based on first-order perturbation 
theory.'*"^) It should be noted that this quantity plays a crucial role in the extraction 
of the second virial coefficient from experimental light-scattering data, due to the 
necessity of extrapolating to zero concentration [141]. 

The "hit-or-miss" algorithm (Section 3.5) can easily be generalized to compute 
the n^^ virial coefficient, given n independent SAW samples. Moreover, the effi- 
ciency should remain reasonably good whenever hyperscaling holds. We do not 
know whether the Barrett and Karp-Luby algorithms can be generalized to virial 
coefficients of order n > 3. Finally, the deterministic Fourier method (Section 3.3) 
can be applied to the third virial coefficient and to at least some of the graphs for 
virial coefficients of order n > 4 (namely, those graphs that can be decomposed 
into series and parallel connections). For example, to compute B3 one first uses the 
Fourier method to compute pi2{x) and thus /i2(x), and likewise for 13 and 23; one 
then computes (/12 * /13 * /23)(0) by passing to Fourier space. 

Another interesting extension of our work would be to study the Domb- Joyce 
model [177] of weakly self-repelling walks, in which each self-intersection is penalized 
by a weight (0 < A < -|-oo); obviously this model interpolates between ordinary 
random walks (A = 0) and SAWs (A = -|-oo). This model can be studied in two 
very different limits: 

(la) iV ^ 00 at fixed A > 0. 

(II) iV ^ 00, A ^ with X = AiV2-'i/2 fixed. 

Case la is the good-solvent regime, and the universal quantities should take exactly 
the same values as for SAWs (for any A > 0). On the other hand, the non-universal 

■*^See DesCloizeaux-Jannink [12], p. 544, equation (13.1.28). Note that in their notation B3/B2 = 
^h{z)/g{zY [compare their equation (13.1.2)]. 
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quantities will manifestly be A-dependent; we expect the approach to to be from 
below for small A, and from above for large A. In particular, there will be some 
intermediate value A* for which the leading correction to scaling will vanish. (In 
Figure 10 this would be achieved e.g. by a family of Hamiltonians that crosses the 
critical surface about halfway between P and Q.) Of course, subleading corrections 
to scaling will still be present, but these will be suppressed by iV"^^, which decays 
much more rapidly than N~^^ . By studying the Domb- Joyce model systematically 
in the (A, iV)-plane, it should be possible to locate A* empirically, and to exploit this 
knowledge to obtain estimates of i^, \E'* and other universal quantities that are less 
contaminated by the effects of corrections to scaling. (We thank Jim Barrett and 
Bernie Nickel for this observation [92].) This is in the same spirit as the analyses of 
series expansions by Chen, Fisher and Nickel [38,41], who found that by varying an 
irrelevant parameter and imposing universality, they could reduce the sensitivity of 
exponent estimates to the effects of corrections to scaling. 

Case II in the Domb- Joyce model should be given (rather trivially) by the con- 
tinuum Edwards model. This fact is hardly in doubt, but the results could serve 
as a useful nonperturbative check on the reliability of the extrapolations to large z 
in the Edwards model [90,91]. This is a warm-up for the problem of the crossover 
scaling behavior near the theta point. 

The deepest — and most difficult — extension of this work would be to SAWs 
with nearest-neighbor attraction, for the purpose of studying the crossover scaling 
behavior near the theta point. In dimension c/ = 3, the crossover scaling functions 
are predicted [14,15] to be given exactly, modulo two nonuniversal scale factors, by 
the continuum Edwards model. In dimension c/ = 2, there is as yet no theoretical 
prediction for the crossover scaling functions, but there are some predictions for 
critical exponents [178]. It is a highly nontrivial problem to develop Monte Carlo 
algorithms that work well near the theta point. The pivot algorithm (taken alone) 
does not work terribly well in this regime [179,180]. 

A Some Geometrical Theorems 

In this appendix we prove some geometrical bounds on T(A, B) = ^(A — B) and 
U(A^B) [see (3.1) and (3.9)] that played a role in Sections 2 and 3 of this paper. 
When averaged over pairs of independent SAWs of length iV, these bounds show 
that 

(a) (T) < const X N^" 

(b) (T) > const X N^" 

(c) {U) > 21ogiV 

modulo the usual assumptions on the scaling of individual SAWs. In particular, 
bounds (a) and (b) together prove (modulo these assumptions) hyperscaling for 
SAWs in dimension d = 2. 
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A.l Theorems and Proofs 

Let us begin by defining several measures of the size of a set A C 1'^: 

• The span in the a^^ coordinate direction 

Sa{A) = maxe„ • x — mine„ • x , (^-1) 

where e„ is the unit vector along the a axis (!<«<(/). 

• The Euclidean diameter 

diam(A) = max \\x — x'\\2 , (^-2) 

x,x'eA 

where || • ||2 denotes the Euclidean norm. 

• The sup-norm diameter 

diamoo(A) = max \\x — x'\\co = max Sa{A) , (^-3) 

where 11 

l<a<d 

Let i? be a positive real number. A set A C Z'' is said to be R-connected if for 
all ' G A there exists ' of points in A such 

that llxj — Xi-i ||oo < -R for 1 < z < n. A set is said to be connected if it is 1-connected. 
(Note that this allows diagonal as well as nearest-neighbor connections.) 

Finally, we let Ad(-R) be the maximum number of lattice points in a closed 
Euclidean ball of radius i?, namely 

Ad(i?) = sup G Z'^: - II2 < i?} • (A.4) 

Note that Ad(i?) = ts,R^ + 0{R^-^) as i? ^ 00, where = 27r'^/V[c/r(c//2)] is the 
volume of the unit ball in R''. 

Theorem A.l Let A,B dT^ . Then 

#{A-B) < n + + 1] (A.5) 

a = l 

#{A-B) < Ad(diam(A) + diam(5)) (A. 6) 



Theorem A. 2 Let A^B he R-connec ted subsets of (d>2). Let Gd C 0(d) be 
the orthogonal symmetry group of~L.^ [so that ^{Gd) = 2'^d\J. Then 

— ^ ^ #{A-gB) > -^dmm^{A)dmm^{B) . (A.7) 
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Theorem A.3 Let A,B C . Then 

J^m,n-l 1 Xf _ . _|_ 1 

U{A,B) > 2 J2 T + . ; (A.8a) 

]^ min 

> log{Af^a. + 1) + log AC.™ , (A.8b) 

where Mmin = ^(B)] and Mmax = ni^'^iT^l^), t^(-^)]- /ir^i inequality 
is the best possible in terms of Mmin o-i^d Mmax- 



Proof of Theorem (A.l). These bounds follow immediately from the fact that 
the set A — B can be contained in a rectangular parallelopiped or a sphere of the 
indicated size. ■ 



Proof of Theorem (A. 2). Let us consider the case d = 2 first. Fix g G G2. 
There exist / > 1 and Xi, . . . , G A such that \\xi — Xi^i | |oo < -R for i = 2, . . . , / 
and I |x/ — Xi I |oo = diamoo(A). Similarly, there exist J > 1 and j/i, . . . , yj G gB such 
that IIj/j - yj-i||oo < -R for j = 2, . . . , J and \ \yj - yi||oo = diamoo(y-B) = diamoo(-B). 
Then Lemma 8.1 of Aizenman [111] implies that 

#{A-gB) > llxj - xilla ||yj - yi||2 I sin6'(a;/ - a;i,yj - yi)| 

> — diamoo( A) diamoo(-B) I sin 6'(a;/ - xi, yj - yi) I (A. 9) 

where 6{v,w) is the angle that the vector w makes with respect to the vector v. 
Let y_|_ G G2 be the operator of rotation by +7r/2. Then for any vectors v and 
6{v, g^w) = 6{v, t/;) + 7r/2, so 

\sin6{v,w)\ + I sin ^(u, y_|_t/;)| = \sin6{v,w)\ + \cos6{v,w)\ > 1. (A. 10) 

From (A. 9), we see that 

#(A - gB) + #(A - g+gB) > diamoo(A) diamoo(5) . (A.ll) 

The theorem for d = 2 now follows by averaging (A.ll) over all y G G2 and dividing 
both sides by 2. 

Now consider the case d > 3. Fix g G Gd- Choose indices /3 and /3' so that 
Sf3(A) = diamoo(A) and Sf)i{gB) = diamoo(-B). If /3 7^ /?', then let P be the 
orthogonal projection of Z onto the (/3,/3') coordinate plane. If /3 = then let 
fi" be any element of {1, . . . , c/} \ and let P be the orthogonal projection onto 
the (/3,/3") coordinate plane. Let h G Gd be the rotation by +7r/2 in the {fi^fi') [or 
(/3,/3")] coordinate plane, which leaves all the other coordinates fixed. Note that 
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Ph = g+P^ where G G2 was defined above. We now use the resuh (A. 11) for 
d = 2: 

#(A - gB) + #(A - hgB) > #{PA - PgB) + #(PA - PhgB) 

> — diamoo(-P-4) diamoo(-Pfi'-B) 
= 4?diamoo(A) diamoo(5) . (A. 12) 

The theorem now foUows by averaging, as for the d = 2 case. ■ 

We remark that some strengthened versions of Theorem A. 2, while quite plausi- 
ble at first sight, are wrong. For example, we at first thought that for 1-connected 
sets (and in particular walks) in Z^, 

#{A-B) > const X Si{A)S2{B). (A.13) 

This is false, since if A and B are both iV-step self-avoiding walks from (0, 0) to 
(iV/2,iV/2) that stay in the diagonal strip \xi — < 1, then ^{A — B) ~ N while 
Si{A)S2{A) ~ iV^. Thus one must average over rotations in some way. We had also 
conjectured that for any sets A,B C. Z^, 

#(A -B) + #iA - g^B) > C #Proj,(A) #Proi,iB) , (A.14) 

where Proj^ is the projection onto the Xa coordinate axis, and g^ is rotation by 7r/2; 
here C is a constant that is independent of the connectedness R. A counterexample 
to this conjecture was found for us by G.L. O'Brien, as follows: Let M be a large 
integer, and let 

A = B = {(0, 0) + k{l, M) + /(M, -1) : < A;, / < M} . (A. 15) 

(This is a large sparse square, slightly tilted.) Then ^Proj;^(A) = ^Proj2(-B) = M^, 
but #(A - B) = #(A - g+B) ^ AM^. We stiU do not know if (A.14) holds if its 
left-hand side is replaced by the left-hand side of (A. 7) with d = 2. 

Proof of Theorem (A. 3). Throughout this proof, we shall assume without loss 
of generahty that M = #(A) > = #{B). 

To show that the first inequality is the best possible, let A (resp. B) consist J\f\ 
(resp. A/2) consecutive sites on the Xi axis (so A and B are parallel rods). Then it 
is easy to check that the first inequality holds as an equality. 

Let u be a vector in R such that the inner-product mapping x ^ v ■ x \s 
one-to-one on AU B. (For example, the coordinates of v could be c/ irrational 
numbers that are linearly independent over the rationals.) Denote the elements of 
A by ai, . . . , a^^, ordered so that v ■ 04 < v ■ a^+i for every i. Similarly, denote the 
elements of i? by 61 , ... , hj^^, ordered so that v ■ < v ■ bj^i for every j . 

We shall use the following observation: 

Observation 1 : If — bj^^ = ai — for some z, /, m, then I < i. 
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(To prove this, note that the hypothesis imphes that — ai = bj^^ — 6^, and so 
V ■ ai — V ■ ai = V ■ bj^^ — v ■ b^ > 0. Therefore v ■ ai > v ■ ai, which imphes / < i.) The 
proof of the foUowing observation is exactly analogous: 

Observation 2 : If aj^^ — bj = — bi for some j, /, m, then I < j ■ 

We now claim that for any z, 

p{a^ -b^r^) < I. (A. 16) 

To prove this, observe that p{x) equals the number of distinct a;'s in A such that 
X G ai — B. But if — bj^^ G ai — then / < z by Observation 1. This proves 
(A. 16). Next, we can use Observation 2 in an analogous way to prove that for any 
J, 

p{aM, - bj) < J. (A.17) 
Now define the following subsets of A — B: 

S' = {a,-bM,: I = l,...,Af2-l} (A.lSa) 
S" = {a^^-bf. J = l,...,M2-l} (A.lSb) 
S'" = {afc-6AA,: A; = Ar2,...,M} (A.lSc) 

It follows from Observation 1 that S' and S" are disjoint, and from Observation 2 
that S" and S'" are disjoint; trivially S' and S'" are disjoint. Therefore 

where the last inequality uses (A. 16), (A.17), and the trivial bound p(x) < M2. A 

n 

simple comparison of Riemann sum to integral shows that > log(n + 1). ■ 

k=i 



A. 2 Application to SAWs 

Now let us apply these bounds to the case in which A and B are independent 
SAWs u^^\u^'^'^ of length N . We first note that that the spans Si, . . . , Sd and the 
diameters D and D^o of either one of these SAWs are expected to scale like iV in 
the sense that 

for any exponents ki, . . . , A;^,^, m > 0. (In particular, this holds if the probability 
distribution of the SAW, with lengths rescaled by iV, converges weakly to some 
probability measure on a space of continuum chains with respect to which the spans 
have finite nonzero moments.) 
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Assuming such scaling, it follows that \I/jv is bounded as iV — )• oo. PROOF: By 
Theorem A.l we have 

T{J'\J'^) < li + + ; (A.21) 

a = l 

but from this and (A. 20) it easily follows that 

^ = {T{J^\J^^))n,n < const X iV'^'' . (A.22) 

Q.E.D. Assuming the usual scaling (2.6) for the radius of gyration, it follows that 
\I/jv is bounded as iV — )• oo. 

Equation (A. 20) also implies that 

{T{J^\J'^^))n,n > const xiV^'' (A. 23) 

in any dimension d > 2. This may be deduced from Theorem A. 2 as follows. Since 
Sn (the set of all iV-step SAWs which start at the origin) is invariant under lattice 
symmetries, we have 

(T(cu(^),cu(2)))^.,^. = ^ ^ #{J'^-J'y) 



> -^j- X] diamoo(cu(^)) diamoo(cu(^)) 
1 



2r2 



3(^00)^, (A.24) 

and combining this with (A. 20) gives (A. 23). 

Lastly, we observe that the scaling assumption (A. 20) implies that hyperscaling 
holds in two dimensions. Indeed, equations (A.22) and (A. 23) together imply that 

{T{ J^\ J^'>))n,n ^ const X N^"" when c/ = 2. (A. 25) 



B Adequacy of Thermalization in the Pivot Al- 
gorithm 

As pointed out in Section 2.2, the initialization of the pivot algorithm is a highly 
nontrivial issue. For not-too-large N (up to a few thousand at least), one can use the 
dimerization algorithm [119,1,81] to produce a perfect equilibrium start. However, 
for very large iV, dimerization is unfeasible, and it is necessary to "thermalize" the 
system by discarding the first ridisc ^ Texp ~ N/f iterations. But this is painful, 
because Texp is a factor ~iV larger than Tint,A for global observables A, and for very 
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large N (> 10^) the CPU time of the algorithm could end up being dominated by 
the thermalization. One is therefore tempted to cut corners in the choice of ridisc- 
In this appendix we want to illustrate how acquiescing to this temptation can lead 
to disaster in the form of large systematic errors. 

Figure 11 shows the temporal history of our pivot-algorithm run for the 3- 
dimensional SAW at iV = 70000. The initial configuration was a pair of parallel 
rods. We averaged the observables and T over bins of 10^ pivot iterations each, 
and plotted the resulting averages for the first 40 bins. Clearly, severe initialization 
bias is present in Rl until at least time 4 X 10^, and in T until at least time 8 X 10^. 
(Here we used the following seat-of-the-pants criterion: if the value at bin n is 
larger than all subsequent bins and smaller than all preceding bins, then severe 
initialization bias is present at time n.) Significant initialization bias — enough to 
cause a systematic error comparable to the (very small) statistical error — could 
well be present at times twice or three times this. So, to be safe, we discarded in 
this case the first 3 X 10^ iterations. (The total run length was approximately 8 X 10^ 
iterations, so a little less than half of the run was discarded.) 

In general, we expect the thermalization to require a time proportional to Texp ~ 
iV//, where / is the pivot-algorithm acceptance fraction. From the run shown 
in Figure 11, we infer that severe initialization bias is present until at least time 
~ 3N/f. Therefore, to be safe, we have always discarded at least 9N/f iterations 
(except of course for runs using a dimerized start). We believe that this rule is 
sufficiently conservative to render the systematic errors arising from inadequate 
thermalization much smaller than the statistical errors. 



C Some Statistical Subtleties 

Let Ai, . . . ^ An be the time series for some observable A (in equilibrium). As 
discussed in [10, Appendix C], the error bar on the sample mean A = ZlILi ^« 
is 

stddev(A) ^ (^ ^^^n^^ACAAiO) ^ ^ 

where 

CAAit) = {A,A,+,) - {Ay (C.2) 

PAAit) = CAAit)/CAAiO) (C.3) 

-j^ CO 

Tmt,A = PAA{t) (C.4) 
^ t=-oo 

Since Tint,A and Caa{^) are not known, they must be estimated from the time series: 
we define 

CAA(t) ^ ^^j:(A.-A)(A+\t\-A) (C.5) 

I I 8 = 1 
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^p^Jt) ^ CAA{t)/CAA{0) (C.6) 
I M ^ 

T,nt,A = J2 PAAit) (C.7) 

^ t=-M 

where we must still choose the "window" M . 

One reasonable approach to choosing M is the "automatic windowing algorithm" 
of [10, Appendix C]: choose M to be the smallest integer such that M > CTint,A{M)^ 
for a suitable "window factor" c. If pAA{t) were roughly a pure exponential, then it 
would suffice to take e.g. c 5 (since < 1%). However, for global observables 
A in the pivot algorithm, pAA{t) is expected to be very slowly decaying: after a 
brief initial decay, one expects pAA{t) ~ up until a time of order Texp ~ N/ 

after which time pAA{t) decays rapidly. This is the behavior exhibited by the ex- 
act solution of the pivot algorithm for ordinary random walk, where pAA{t) ~ 1/t 
in the intermediate region [10, Section 3.3]. We confirmed this behavior empir- 
ically for the self-avoiding walk by making an extremely long (n = 10® pivots) 
simulation at iV = 1000 (d = 2): the sample autocorrelation functions Paa{^) 
A = Rl^ R'^^TKarp-Luby,R=20 s.Te showu in a log-log plot in Figure 12. For there 
is a wide intermediate region (20 ^ t < 1000) in which the log-log plot is roughly 
straight, yielding pAA{t) ~ 1/^' with q 1.0-1.2. For there is more curvature, 
but the slope q is in the same ballpark. For TKarp-Luby,R=20^ in contrast, the curva- 
ture is so great that we are unable to identify an intermediate region with clear 1/t' 
behavior. In any case, we have rough empirical confirmation for the theoretically 
predicted behavior; we are unable to say whether q = lovq = l-\-pov something 
else, but clearly g 1. 

Because of this very slow decay, the automatic windowing algorithm leads to 
significant underestimates of Tint,A^ even when the window factor c is as large as 
10 or 20. We therefore defined a modified estimator Tint,A by extrapolating pAA{t) 
proportionally to 1/t in the region beyond the window, i.e. Tint,A = Mpj^j^(M)/t for 
t > M, and cutting off the sum at a much larger time M* of order Texp- 

-r M 

TraiA = H PAA(t) + M^p ^JM) log(M* / M) . (C.8) 

^ t=-M 

(Since M^M* ^ 1, we approximated the second sum by an integral.) Here M is 
defined, as before, by the automatic windowing algorithm. In the absence of any 
precise knowledge of Texp^ we took M* = N/f; luckily Tint,A is not very sensitive to 
the choice of M*, because of the logarithmic dependence and because /9^^(M) <^ 1. 

In Table 10 we show the results for a typical run (c/ = 2, iV = 80000) for the 
observables A = Rl, R^, Rl^,TKarp-L uby,R=i50 as a fuuctiou of the window factor c. 
Time is measured in units of 25 pivots. Note that: 

(a) The standard windowing estimates f 

int. A 3, factor 2-3 smaller than the 
modified estimates Tint,A7 even at c = 20-30. 
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(b) The estimates Tint,A are slowly but relentlessly increasing as a function of c; 
it is perfectly plausible that they will double or triple by the time c reaches 
infinity. 

(c) The estimates Tint,A are roughly stable as a function of c as soon as c > 4. 
However, they do show some fluctuation, because they are based on extrapo- 
lation from a single noisy point /9^^(M). 

We therefore chose as our flnal estimate the average of Tint,A for the eleven values 
c = 20, 21, ... , 30; by taking the average, we reduce the fluctuations mentioned in 
(c). 

This whole procedure is, of course, inelegant and ad hoc. But it does work 
reasonably well: we expect that the estimates of Tint,A are accurate to about 10%. 
This is not good enough for a serious study of dynamic critical behavior; but it is 
good enough for our present purpose, which is merely to set error bars on the static 
quantities A. In the future we hope to devise better methods for analyzing time 
series with slow decay of the autocorrelation function. 

D Remarks on the Field-Theoretic Estimates of 
Universal Amplitude Ratios 

The critical exponents and universal amplitude ratios associated with polymers 
in a good solvent can be extracted from any family of theories which intersects 
transversally the domain of attraction (= stable manifold) of the good-solvent flxed 
point Hq^. One such family — which has no special status at Hq^^ but is computa- 
tionally convenient — is the continuum Edwards model (or what is equivalent, the 
n = continuum (p"^ fleld theory) [12]. This model becomes critical (i.e. crosses the 
stable manifold of Hq^) when its bare coupling constant (self-avoidance parameter) 
z tends to oo. 

There are two main approaches to the quantitative study of the continuum 
Edwards model in the limit z ^ oo: 

• Perturbation expansion in the coupling constant z at fixed dimension d = 3 
(or d = 2). In this case the problem is to estimate the asymptotic behavior 
as 2; — )• oo from the flrst few terms of a perturbation expansion around z = Q. 

• Expansion in e = A — d. Here the critical exponents and limiting amplitude 
ratios, which correspond to the limit z — )• 00, can be obtained directly from a 
suitable renormalization-group analysis. The problem is then to extrapolate 
to e = 1 (or e = 2). 

In this appendix we want to summarize the results for the universal amplitude ratios 
K = 6(-R^)/(-Rg), \E'* and \E'^ which have been obtained by each of these methods, 
and to make some comments on their extrapolation. 
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1) Perturbation expansion at fixed dimension d = 3. Let aR(z)^ as{z) and h(z) 
be the conventional expansion and second virial factors of the continuum Edwards 
modeL One then has 

Hz) = aUz)/aUz) (D.l) 
^{z) = zh{z)/ai{z) (D.2) 
^n{z) = zh{z)/aUz) (D.3) 

The known perturbation series in d = 3 foUows [89,107,136,181,91]: 

2 / 32 977r 



^z) = 1 z + 



2 



Z 



35 V105 1296. 
+ 0.056588^3 _ 12202^^ + 0{z^) (D.4) 

^(z) = z — — ^^^^2 ^ 25.5896415457^3 ^ Q/^4^ 

105 

<i>R(z) = z - 254- 128a/2 ^2 ^ 26.1049923667^3 ^ q/^a p 

15 

Crude estimates of = K(oo) can be obtained from the [1/1] and [2/2] Pade 
approximants 

1 _|_ ( 554 _ 33957r \ 

= (D-7) 

L^'-^J^ ^ -, , / 16 33957r^ ^ ^ 



1 _|_ ib _ 33957r 
V 3 2592 J 

1 + 8.08283^ + 10.5132^2 

^f2/2i(^) = ^ ^ ^ (D-8) 

^ ' ^ 1 + 8.13998^ + 10.9087^2 v ; 

yielding ^[i/i] ~ 0.9531 and ^[2/2] ~ 0.9637. More sophisticated extrapolations have 
been performed by Shanes and Nickel [136], who find = 0.9607 ± 0.0018. (Our 
Monte Carlo value is 0.9594 ± 0.0012.) 

Likewise, crude estimates of = \I/(oo) and = "^^(00) can be obtained 
from the [1/1] Pade approximant 



s 



L^'-^J^ ^ -r _|_ 1769-896v^ ^ ^ ^ 

^ ^ 105 ^ 



^^r^,^,(^) = — (D.IO) 

"-'i'-l '-i'- > -, , 254-128v^ ^ ^ ^ 

15 ^ 

yielding ~ 0.2092 and ^'^,[1/1] 0.2055. (Our Monte Carlo values are 0.2471 + 

0.0003 and 0.2322 + 0.0003, respectively.) 

We can also try the "direct renormalization" approach of des Cloizeaux, Conte 
and Jannink [90]: define the "effective exponent" or{z) = (z /2) dlog aj^(z)/ dz^ 
which approaches the limiting value a'^ = 2iy — 1 as 2; — )• 00; re-express \I/ and 
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as functions of aji; and extrapolate these series to aji = cr|j, using the best estimate 
of i^. If we carry out this last step by the most naive method imaginable — namely, 
straight evaluation of the cubic polynomial at cr|j = 2(0.5877) — 1 = 0.1754 — we 
obtain = 0.2483 and = 0.2379, which is not bad at all for such a short series. 

2) Expansion in e = A — d. The limiting universal ratios K*, and can be 
evaluated in dimension d = A — ein powers of e; the results are'*'*: 

= 1 - — - 0.0306286^ + Oie^) (D.ll) 
96 ^ ^ ' ^ 

= - + —(— + \og2\ + Oie") (D.12) 
8 16 V48 ^ ; ^ ^ ^ ' 

= 8 + T^b + ^°S2) (D.13) 

Evaluating these at e = 1, 2, 3 yields 

^ 0.9590, 0.8567, 0.6931 (D.14) 
^* ^ 0.2686, 0.8243, 1.6672 (D.15) 
^ 0.2660, 0.8139, 1.6438 (D.16) 

(We remark that the [1/1] Pade approximants for all three quantities have poles 
between e = and e = 1, and so are unreliable for e > 1.) 

Des Cloizeaux [182, footnote 7]'*^ has suggested to augment the e-expansions 
(D.11)-(D.13) by enforcing the known exact values at d = 1 (e = 3), which are 
K = 1/2, ^* = 2/7r*/2 ^ 1.1284 [see (2.20)] and = (2/7r)*/2 ^ 0.7979. This 
produces the cubic polynomials 

= 1 - — - 0.0306286^ - 0.0071526^ (D.17) 

96 ^ ^ 

e e2 /77 , ,\ / 2 327 9 



= - + — — + log2 + — — log 2 (D.18) 

8 16 V48 ^ y 27 V7ri/2 256 16 ^ ' ^ ^ 

e A25 , \ , 321 9 , \ , 

n = - + - - + log2 + - (2/7r)*/2- — --log2 (D.19) 



16 V16 ° J 27 V ' ' 256 16 
Evaluating this at e = 1,2 yields 

^ 0.9518, 0.7994 (D.20) 

^* ^ 0.2486, 0.6647 (D.21) 

^ 0.2346, 0.5633 (D.22) 



■^■^The expansion for is from [153]; see also equation (13.1.60) of [12]. The expansion for is 
from equation (5.35) of [182]; see also equations (12.3.102)/(13.1.11) of [12]. (In these works, 
is called g* . To establish the connection, compare equations (5.1.107) and (13.1.2) of [12].) The 
expansion for ^* can then be derived from these two. Alternatively, it can be found as the limit 
rj ^ CO of equation (4.15a) of [183]. 

45See also [12, pp. 541 and 557]. 
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For this modification has actuaUy worsened the agreement with the Monte Carlo 
value, both in d = 3 and in d = 2. On the other hand, for and\I/^ this "modified 
e-expansion" prediction is amazingly close to the correct value, both in d = 3 and 
d = 2. It would be useful to obtain a better understanding of whether this is a 
coincidence or not — perhaps by calculating the 0(e"^) term in 
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N 


(T> 


(T2)/104 


var(T)/104 


(U) 


(VBarrett)/10'* 


(Vifarp-Ltiij/)/10'* 


100 


1373.53 (0.48) 


193.95 


5.29 


339.07 


157.22 


483.94 ( 0.28) 


200 


3854.69 (1.39) 


1529.10 


43.23 


723.10 


1435.55 


4236.24 ( 2.58) 


300 


7065.32 (2.68) 


5139.55 


147.68 


1124.57 


5196.87 


15025.26 ( 9.79) 


400 


10860.24 (4.11) 


12146.26 


351.79 


1535.59 


12898.02 


36803.23 ( 24.10) 


500 


15152.67 (8.32) 


23639.57 


679.22 


1950.72 


26002.92 


73527.83 ( 70.29) 



Table 1: Quantities relevant to the Barrett and Karp-Luby algorithms, as a function 
of Ni = N2 = iV, for 2-dimensional self-avoiding walks. Error bars are one standard 
deviation. 



N 


(T) 


(T2)/104 


var(T)/104 


(U) 


(VBarrett)/10'* 


(Vifarp-Ltiij/)/10'* 


100 


3346.95 ( 0.93) 


1136.81 ( 0.63) 


16.60 ( 0.88) 


1719.60 ( 0.75) 


617.35 ( 0.15) 


2337.71 ( 0.90) 


200 


11159.79 ( 7.25) 


12646.28 ( 16.27) 


192.20 ( 22.94) 


5149.61 ( 5.25) 


8158.67 ( 5.25) 


28739.30 ( 26.18) 


500 


55272.12 (17.35) 


310214.61 ( 193.40) 


4713.83 ( 272.40) 


22089.93 (11.11) 


244244.80 ( 87.33) 


795935.08 ( 364.63) 


800 


125986.30 (40.59) 


1611900.23 (1029.67) 


24645.50 (1451.35) 


46787.65 (24.22) 


1389999.88 ( 532.61) 


4381808.49 (2079.36) 


1000 


186508.58 (61.41) 


3532249.52 (2309.75) 


53704.50 (3253.13) 


66896.98 (35.56) 


3170834.30 (1271.38) 


9855248.60 (4819.17) 



Table 2: Quantities relevant to the Barrett and Karp-Luby algorithms, as a function 
of Ni = N2 = iV, for 3-dimensional self-avoiding walks. Error bars are one standard 
deviation. 
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N 

Too" 

200 
300 
400 
500 
1000 
1500 
2000 
2500 
3000 
3500 
4000 
4500 
5000 
5500 
6000 
6500 
7000 
7200 
7500 
8000 
8250 
8400 
8650 
9000 
9200 
9500 
9650 
9700 
9750 
9850 
9900 
10000 
15000 
20000 
40000 
60000 
80000 



777. 
2194. 
4017. 
6185. 
8621. 
24459. 
44774. 
68929. 
96755. 
127104, 
159980, 
194741, 
232822, 
272893, 
315118, 
359467, 
404321, 
450254, 
471499, 
500172, 
551978, 
577862, 
593186, 
620944, 
657926, 
677758, 
710439, 
728607, 
736619, 
738974, 
755190, 
760424, 
770061, 
1418675, 
2177520, 
6177080, 
11352100, 
17390000, 



690 ± 
000 ± 
710 ± 
440 ± 
527 ± 
300 ± 
300 ± 
400 ± 
100 ± 
073 ± 
000 ± 
264 ± 
000 ± 
000 ± 
000 ± 
102 ± 
000 ± 
000 ± 
566 ± 
000 ± 
094 ± 
000 ± 
000 ± 
000 ± 
882 ± 
000 ± 
000 ± 
000 ± 
000 ± 
000 ± 
000 ± 
333 ± 
308 ± 
444 ± 
000 ± 



0.757 
2.180 
4.179 
6.388 
10.642 
42.838 
85.796 
131.083 
216.912 
176.229 
390.518 
327.675 
662.439 
709.164 
821.224 
613.861 
1028.130 
1217.260 
762.481 
1396.910 
926.915 
1508.860 
1532.640 
1835.070 
1129.894 
1865.980 
2295.450 
1967.000 
2076.780 
1968.260 
2200.680 
1121.483 
1294.861 
1916.596 
5305.390 
000 ±17206.500 
000 ±31723.800 
000 ±42700.700 



109, 
307, 
563, 
867, 
1210, 
3425, 
6282, 
9653, 
13559, 
17802, 
22420, 
27332, 
32693, 
38287, 
44140, 
50437, 
56654, 
63233, 
66130, 
70329, 
77426, 
80968, 
83241, 
87038, 
92273, 
95064, 
99690, 
102546 
103406 
103711 
106059 
106678 
108002 
198953 
305098 
864495 
1595490 
2438490 



142 ± 
451 ± 
694 ± 
673 ± 
549 ± 
590 ± 
270 ± 
390 ± 
200 ± 
697 ± 
700 ± 
984 ± 
300 ± 
900 ± 
800 ± 
817± 
600 ± 
900 ± 
359 ± 
000 ± 
100 ± 
200 ± 
100 ± 
900 ± 
834 ± 
200 ± 
100 ± 
000 ± 
000 ± 
000 ± 
000 ± 
709 ± 
619± 
381 ± 
000 ± 



0.060 
0.174 
0.334 
0.512 
0.948 
5.623 
11.785 
18.207 
30.191 
23.913 
52.365 
46.167 
84.724 
106.209 
116.423 
84.619 
149.773 
187.088 
104.382 
207.664 
127.431 
206.938 
213.843 
249.425 
161.439 
270.506 
316.909 
253.146 
304.886 
281.340 
272.353 
159.471 
175.686 
265.417 
771.573 
000 ±2482.060 
000 ±4517.620 
000 ±6045.320 



(Rm) 



10742.500 ± 
19684.600 ± 
30289.000 ± 



207146 
220086 
242583 
253902 
260645 
272775. 
289167 
298003 
312471 
320619 
324039 
325021 
331942 
334263 
338560 
624522 
957108 
2712030 
4991160 
7642600 



970 ± 
000 ± 
184 ± 
000 ± 
000 ± 
000 ± 
124 ± 
000 ± 
000 ± 
000 ± 
000 ± 
000 ± 
000 ± 
606 ± 
337 ± 
711± 
000 ± 
000 ± 
000 ± 
000 ± 



17.119 
34.762 
53.735 



55851.100 ± 88.597 

85730.297 ± 178.520 

119932.000 ± 303.951 

158072.674 ± 255.803 

177615.000 ± 443.270 



379.611 
587.985 
384.604 
602.582 
621.496 
753.115 
464.199 
756.903 
877.201 
781.385 
859.093 
810.977 
854.863 
467.136 
522.242 
1018.853 
2249.350 
7031.920 
13044.200 
17427.300 



1373, 
3854, 
7065, 
10860, 
15154, 
42853, 
78573, 
120773 
169394 
222357 
280117 
341511 
408602 
478110 
551118 
630667 
708214 
790171 
825946 
878189 
967214 
1012461 
1040352 
1087701 
1152906 
1189077 
1247906 
1279259 
1293086 
1297868 
1322669 
1332241 
1350156 
2486038 
3817040 
10807802 
19893345 
30476824 



53 ± 
69 ± 
32 ± 
24 ± 

26 ± 

78 ± 
80 ± 

27 ± 
08 ± 

79 ± 
27 ± 
22 ± 

67 ± 

68 ± 



0.48 
1.39 
2.68 
4.11 
7.85" 
57.45 
115.01 
182.27 
303.43" 
246.92" 
555.52" 
447.94" 
848.48" 
996.75" 
77 ±1220.07" 
93 ± 852.57" 
98 ±1548.77" 
67 ±1863.89" 

29 ±1038.17" 
91 ±2013.06" 
84 ±1288.57" 
00±2061.87" 
40±2145.18" 

30 ±2552.58" 
64 ±1583.77" 
20 ±2590.80" 
30 ±3197.34" 
90 ±2697.09" 
50±2902.13" 
30 ±2753.30" 
90 ±2677.55" 
41 ±1605.57" 
71 ±1853.36" 
05 ± 2606.19 
60 ± 7883.18 
00 ±23097.90 
00 ±48018.50 
00±58801.80 



0.25378 ±0.00024 
0.23452 ±0.00024 
0.22140±0. 00023 

0.20509±0. 00019 

0.19376±0. 00023 



0.17953±0. 00021 

0.17309±0. 00018 
0.16939±0. 00017 

0.16559±0. 00017 



0.16247±0. 00016 
0.15029±0. 00011 
0.14216±0. 00019 
0.12451 ±0.00019 
0.11573±0. 00016 
0.10912±0. 00013 



Table 3: The results of our runs in dimension d = 2. Errors are it one standard 
deviation. " indicates a possible minor bug in measurement; see text. 
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N 








100 





14034 ±0 


00021 


— 





66764 ± 0.00060 


200 





14013 ±0 


00022 


— 





66514 ± 0.00062 


300 





14030 ±0 


00023 


— 





66495 ± 0.00065 


400 





14028 ±0 


00023 


— 





66402 ± 0.00064 


500 





14041 ±0 


00028 


— 





66413 ±0.00086" 


1000 





14005 ±0 


00048 


0.43920 ±0.00147 





66367 ± 0.00198 


1500 





14031 ±0 


00053 


0.43964 ±0.00162 





66353 ± 0.00222 


2000 





14005 ±0 


00053 


0.43942 ±0.00162 





66373 ± 0.00225 


2500 





14014 ±0 


00063 


— 





66277 ±0 


00266" 


3000 





14006 ±0 


00038 


0.43941 ±0.00131 





66262 ±0 


00163" 


3500 





14015 ±0 


00067 


— 





66281 ±0 


00286" 


4000 





14036 ±0 


00047 


0.44023 ±0.00166 





66285 ±0 


00199" 


4500 





14042 ±0 


00076 


— 





66304 ±0 


00310" 


5000 





14030 ±0 


00075 


0.43948 ±0.00226 





66247 ±0 


00322" 


5500 





14008 ±0 


00073 


— 





66237 ±0 


00321" 


6000 





14031 ±0 


00048 


0.43974 ±0.00146 





66335 ±0 


00201" 


6500 





14012 ±0 


00073 


0.43929 ±0.00221 





66318 ±0 


00320" 


7000 





14044 ± 


00080 


— 





66293 ±0 


00353" 


7200 





14026 ±0 


00045 


0.43934 ±0.00152 





66260 ±0 


00188" 


7500 





14061 ±0 


00081 


0.44002 ±0.00240 





66245 ± 


00347" 


8000 





14027 ±0 


00047 


0.43948 ±0.00143 





66273 ±0 


00197" 


8250 





14012 ±0 


00072 


0.43938 ±0.00219 





66338 ±0 


00305" 


8400 





14033 ±0 


00072 


0.43940 ±0.00218 





66304 ±0 


00307" 


8650 





14017 ±0 


00082 


0.43929 ±0.00251 





66297 ±0 


00346" 


9000 





14025 ±0 


00049 


0.43951 ±0.00146 





66285 ±0 


00207" 


9200 





14026 ±0 


00079 


0.43969 ±0.00233 





66358 ±0 


00333" 


9500 





14032 ±0 


00090 


0.43983 ±0.00266 





66409 ± 


00381" 


9650 





14074 ±0 


00073 


0.44004 ±0.00226 





66182 ±0 


00303" 


9700 





14038 ±0 


00081 


0.43990 ±0.00241 





66341 ±0 


00344" 


9750 





14034 ±0 


00075 


0.43983 ±0.00227 





66390 ±0 


00321" 


9850 





14044 ± 


00077 


0.43955 ±0.00241 





66161 ±0 


00304" 


9900 





14029 ±0 


00042 


0.43958 ±0.00126 





66253 ±0 


00179" 


10000 





14025 ±0 


00046 


0.43965 ±0.00142 





66321 ±0 


00199" 


15000 





14024 ±0 


00038 


0.44022 ±0.00131 





66291 ± 0.00158 


20000 





14011 ±0 


00070 


0.43954 ±0.00210 





66372 ± 0.00305 


40000 





13995 ±0 


00079 


0.43905 ±0.00236 





66324 ± 0.00332 


60000 





14055 ±0 


00079 


0.43967 ±0.00238 





66147 ± 0.00347 


80000 





14022 ±0 


00069 


0.43948 ±0.00208 





66305 ± 0.00292 



Table 4: Universal amplitude ratios in dimension d = 2. Errors are ± one stan- 
dard deviation, based on triangle inequality. " indicates a possible minor bug in 
measurement; see text. 
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N 


{Rl) 


(Rl) 


(T) 


/ 


100 


265 


053 ± 





169 


41 


912 ± 





017 


3347 


31 ± 


0.62 





60873 ±0 


00012 


150 


429 


084 ± 





392 


67 


954 ± 





047 


6756 


65 ± 


2.70 





58224 ±0 


00008 


200 


603 


999 ± 





732 


95 


750 ± 





069 


11148 


05 ± 


3.40 









300 


976 


094 ± 





793 


154 


937 ± 





076 


22613 


93 ± 


4.75 









400 


1372 


160 ± 


2 


156 


217 


745 ± 





204 


37357 


58 ± 


14.70 









500 


1786 


826 ± 


1 


266 


283 


943 ± 





131 


55254 


94 ± 


11.95 





51016±0 


00012 


600 


2217 


770 ± 


2 


568 


352 


598 ± 





243 


76076 


21 ± 


21.87 









700 


2661 


450 ± 


5 


081 


423 


525 ± 





676 


99752 


01 ± 


126.77" 









800 


3113 


250 ± 


3 


530 


495 


619 ± 





346 


125986 


30 ± 


40.59 









1000 


4053 


528 ± 


3 


937 


645 


529 ± 





413 


186509 


18 ± 


59.65" 









1200 


5023 


820 ± 


9 


959 


801 


123 ± 


1 


321 


256749 


19 ± 


333.41" 









1500 


6545 


540 ± 


12 


712 


1042 


710 ± 


1 


896 


379923 


40 ± 


529.85 








1800 


8092 


940 ± 


15 


878 


1290 


860 ± 


2 


302 


523305 


41 ± 


726.03" 









2000 


9170 


600 ± 


21 


899 


1461 


760 ± 


3 


060 


629686 


03 ± 


1061.54" 









2500 


11926 


100 ± 


23 


465 


1900 


610 ± 


3 


589 


931639 


59 ± 


1357.98" 









2800 


13632 


500 ± 


29 


044 


2173 


670 ± 


4 


100 


1136711 


90 ± 


1612.15" 









3000 


14791 


272 ± 


16 


145 


2361 


192 ± 


2 


402 


1285631 


38 ± 


1029.93" 





41981 ±0 


00011 


3250 


16289 


800 ± 


28 


272 


2596 


720 ± 


3 


918 


1481692 


70 ± 


1778.01" 









3500 


17762 


000 ± 


37 


294 


2836 


320 ± 


5 


573 


1686417 


60 ± 


2668.40" 









3800 


19526 


200 ± 


41 


537 


3119 


110 ± 


6 


092 


1946853 


70 ± 


3045.89" 









4500 


23857 


422 ± 


38 


095 


3812 


571 ± 


5 


218 


2626477 


19 ± 


2858.05" 









5000 


27006 


559 ± 


29 


538 


4315 


483 ± 


4 


431 


3159180 


21 ± 


2662.07" 





39731 ±0 


00011 


6000 


33467 


400 ± 


43 


910 


5341 


140 ± 


6 


595 


4354063 


40 ± 


4277.04 





38938 ±0 


00010 


7000 


zLfll 9zL 


100 ± 


51 


Q1 Q 




520 ± 


7 


D ( D 




60 ± 


5789.41 








8000 


46951 


672 ± 


57 


753 


7496 


853 ± 


8 


063 


7217286 


50 ± 


6540.09" 








10000 


61128 


600 ± 


88 


087 


9749 


630 ± 


12 


848 


10702551 


00 ± 


12094.20 





36854 ±0 


00011 


13000 


83031 


100 ± 


122 


458 


13264 


800 ± 


19 


003 


16973895 


00 ± 


19439.20 





35820 ±0 


00011 


15000 


98538 


647 ± 


104 


138 


15743 


182 ± 


16 


204 


21865346 


33 ± 


18392.31 





35308 ±0 


00008 


23000 


162265 


000 ± 


254 


103 


25955 


600 ± 


38 


810 


46375397 


00 ± 


59208.20 





33728 ±0 


00011 


25000 


179897 


000 ± 


269 


833 


28727 


500 ± 


41 


726 


53815101 


00 ± 


67927.40 





33454 ±0 


00011 


30000 


222088 


000 ± 


342 


061 


35490 


300 ± 


54 


269 


74087837 


00 ± 


97497.60 





32767 ±0 


00011 


35000 


266714 


000 ± 


432 


848 


42626 


800 ± 


65 


234 


97208448 


00 ± 


119367.00 





32255 ±0 


00010 


40000 


311740 


821 ± 


274 


332 


49801 


999 ± 


41 


934 


122939664 


03 ± 


87490.24 





31793 ±0 


00006 


70000 


602779 


000 ±1304 


950 


96255 


700 ±196 


856 


330049626 


00 ±550259.00 





29983 ±0 


00013 


80000 


704651 


925 ± 


599 


289 


112663 


799 ± 


89 


456 


417196845 


49 ±279956.55 





29561 ±0 


00005 



Table 5: The results of our runs in dimension d = 3. Errors are ± one standard 
deviation. " indicates a possible minor bug in measurement; see text. 
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N 






100 





15813 ±0 


00016 





27693 ± 0.00022 


150 





15837 ±0 


00025 





27077 ± 0.00039 


200 





15853 ±0 


00031 





26710 ± 0.00037 


300 





15873 ±0 


00021 





26323 ± 0.00025 


400 





15869 ±0 


00040 





26100 ± 0.00047 


500 





15891 ±0 


00019 





25924 ± 0.00024 


600 





15899 ±0 


00029 





25794 ± 0.00034 


700 





15913 ±0 


00056 





25691 ±0.00094" 


800 





15920 ±0 


00029 





25632 ± 0.00035 


1000 





15925 ±0 


00026 





25528 ±0.00033" 


1200 





15946 ±0 


00058 





25418 ±0.00096" 


1500 





15930 ±0 


00060 





25330 ± 0.00104 


1800 





15950 ±0 


00060 





25329 ±0.00103" 


2000 





15940 ±0 


00071 





25293 ±0.00122" 


2500 





15937 ±0 


00061 





25240 ±0.00108" 


2800 





15945 ±0 


00064 





25179 ±0.00107" 


3000 





15963 ±0 


00034 





25154 ±0.00059" 


3250 





15941 ±0 


00052 





25137 ±0.00087" 


3500 





15968 ±0 


00065 





25062 ± 0.00114" 


3800 





15974 ±0 


00065 





25088 ±0.00113" 


4500 





15981 ±0 


00047 





25046 ±0.00079" 


5000 





15979 ±0 


00034 





25016 ±0.00060" 


6000 





15959 ±0 


00041 





25040 ± 0.00071 


7000 





15969 ±0 


00040 





24976 ± 0.00070 


8000 





15967 ±0 


00037 





24960 ±0.00063" 


10000 


n 






n 


94.0^7 _i_ n nnn7s 

Zj^vo i nz u.uuuio 


13000 





15976 ±0 


00046 





24941 ± 0.00082 


15000 





15977 ±0 


00033 





24849 ± 0.00059 


23000 





15996 ±0 


00049 





24896 ± 0.00088 


25000 





15969 ±0 


00047 





24811 ± 0.00085 


30000 





15980 ±0 


00049 





24875 ± 0.00090 


35000 





15982 ±0 


00050 





24795 ± 0.00087 


40000 





15975 ±0 


00028 





24832 ± 0.00049 


70000 





15969 ±0 


00067 





24810 ± 0.00117 


80000 





15989 ±0 


00026 





24766 ± 0.00046 



Table 6: Universal amplitude ratios in dimension d = 3. Errors are ± one stan- 
dard deviation, based on triangle inequality. " indicates a possible minor bug in 
measurement; see text. 
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a 


h 


A 


x' 


100 
200 
300 
500 
1000 


0.16006 ± 0.00027 
0.15996 ± 0.00026 
0.15992 ± 0.00024 
0.15985 ± 0.00022 
0.15986 ± 0.00029 


-0.00084 ±0.00032 
-0.00072 ±0.00032 
-0.00067 ±0.00031 
-0.00059 ±0.00030 
-0.00059 ±0.00032 


0.370 ±0.123 
0.466 ±0.227 
0.526 ±0.288 
0.685 ±0.498 
0.668 ±0.950 


2.20 (32 DF, level ^ 100%) 
1.93 (30 DF, level ^ 100%) 
1.78 (29 DF, level ^ 100%) 
1.38 (27 DF, level ^ 100%) 
1.36 (23 DF, level ^ 100%) 



Table 7: Fits {R]) / {Rl) = a + 6(iV/1000)-^ with a, 6, A all variable, for 3- 
dimensional SAWs. The true error bars are probably 1/3 of those indicated 
here (see text). 





^* = a 


h 


A 


x' 


100 
200 
300 
500 
1000 
1500 


0.24710 ± 0.00027 
0.24705 ± 0.00034 
0.24714 ± 0.00036 
0.24717 ± 0.00045 
0.24711 ± 0.00063 
0.24715 ± 0.00083 


0.00819 ± 0.00029 
0.00826 ± 0.00041 
0.00813 ± 0.00046 
0.00808 ± 0.00056 
0.00815 ± 0.00062 
0.00822 ± 0.00089 


0.561 ±0.013 
0.556 ±0.026 
0.569 ±0.035 
0.573 ±0.057 
0.562 ±0.105 
0.574 ±0.201 


3.67 (32 DF, level ^ 100%) 
3.57 (30 DF, level ^ 100%) 
3.22 (29 DF, level ^ 100%) 
3.06 (27 DF, level ^ 100%) 
2.92 (23 DF, level ^ 100%) 
2.82 (21 DF, level ^ 100%) 



Table 8: Fits \I/ = a ± 6(iV/1000) ^ with a,b,A all variable, for 3-dimensional 
SAWs. The true error bars are probably 1/3 of those indicated here (see text). 





d = 1 


d = 2 


(/ = 3 


(/ = 4 


hard— sphere 


1.1284 


4/3 


1.6186 


2 




0.9772 


0.7427 


0.6585 


0.6755 




0.9974 


0.7048 


0.5000 






1.0197 


0.7004 


0.4508 






1.0365 


0.6938 


0.4216 






1.0490 


0.6918 


0.4033 






1.0586 


0.6891 


0.3884 






1.0662 


0.6882 


0.3777 






1.1284 


0.6630 


0.2471 





^ hard— sphere 


1 


0.4972 


0.1527 






Table 9: \I/jv for short chains (rounded to four decimal places) from exact enumer- 
ations [132], along with our Monte Carlo values for 
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A = 


e 


A = 


Rl 






A = 


n2 

m 


A 


= Tk 




c 














Tint,T 


Tint,T 


2 


1 


198 


5.213 


1.873 


11 


247 


1 


455 


8.128 





914 


5.227 


3 


1 


270 


4.667 


2.287 


11 


465 


1 


771 


7.413 


1 


168 


6.564 


4 


1 


363 


4.166 


2.641 


11 


523 


1 


898 


7.153 


1 


350 


7.161 


5 


1 


394 


3.837 


2.889 


11 


489 


1 


997 


7.117 


1 


492 


7.470 


6 


1 


447 


3.914 


3.074 


11 


562 


2 


108 


6.806 


1 


613 


7.905 


7 


1 


487 


3.609 


3.236 


11 


918 


2 


194 


6.547 


1 


713 


7.925 


8 


1 


518 


3.602 


3.374 


12 


221 


2 


241 


6.419 


1 


851 


8.640 


9 


1 


532 


3.526 


3.523 


12 


023 


2 


313 


7.317 


1 


955 


7.801 


10 


1 


555 


3.380 


3.648 


12 


255 


2 


374 


7.144 


2 


051 


8.661 


11 


1 


575 


3.257 


3.755 


11 


394 


2 


433 


7.482 


2 


135 


8.936 


12 


1 


600 


4.291 


3.833 


11 


911 


2 


484 


7.606 


2 


215 


9.243 


13 


1 


612 


4.265 


3.918 


11 


300 


2 


530 


7.276 


2 


280 


8.321 


14 


1 


632 


3.793 


3.999 


12 


654 


2 


571 


7.237 


2 


341 


8.508 


15 


1 


651 


4.066 


4.092 


12 


356 


2 


621 


7.376 


2 


423 


9.482 


16 


1 


672 


4.428 


4.168 


13 


123 


2 


653 


6.712 


2 


476 


8.792 


17 


1 


690 


4.196 


4.232 


11 


731 


2 


684 


7.210 


2 


528 


8.719 


18 


1 


708 


4.113 


4.303 


12 


489 


2 


710 


6.342 


2 


591 


9.608 


19 


1 


723 


4.005 


4.384 


15 


178 


2 


734 


7.033 


2 


653 


9.815 


20 


1 


736 


4.045 


4.446 


13 


261 


2 


770 


7.752 


2 


691 


9.078 


21 


1 


750 


4.531 


4.511 


11 


689 


2 


798 


7.535 


2 


741 


9.663 


22 


1 


763 


4.260 


4.545 


9 


002 


2 


828 


7.803 


2 


785 


7.961 


23 


1 


776 


4.432 


4.587 


10 


681 


2 


853 


7.397 


2 


824 


10.257 


24 


1 


787 


3.839 


4.623 


12 


226 


2 


883 


6.922 


2 


873 


9.043 


25 


1 


799 


4.198 


4.664 


11 


623 


2 


897 


6.365 


2 


917 


10.258 


26 


1 


813 


3.530 


4.708 


12 


612 


2 


913 


7.024 


2 


958 


9.786 


27 


1 


820 


3.389 


4.738 


8 


981 


2 


940 


7.839 


2 


991 


7.904 


28 


1 


827 


4.013 


4.761 


11 


195 


2 


963 


8.450 


3 


024 


9.305 


29 


1 


837 


4.626 


4.783 


8 


318 


2 


993 


8.515 


3 


060 


11.170 


30 


1 


845 


4.474 


4.806 


11 


310 


3 


020 


9.014 


3 


119 


13.602 



Table 10: Standard windowing estimate Tint,A and modified estimate Tint,A for the 
observables A = Rl, R'^, R'^,TKarp-Luby,R=i50, as a function of window factor c, for 
SAWs on the square lattice at iV = 80000. Time is measured in units of 25 pivots. 
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0.752 



0.751 



d = 2 SAW, Pure Power-Law Fit 



0.750 



0.749 



0.748 



10' 




10' 



10' 



Figure 1: Estimated exponent v from pure power-law fits to (-Re)^ {^gY (^) 

and {RmY (o), and estimated exponent (2A4 — 7)/2 from pure power-law fit to (T) 
(*), plotted versus Nmin- Note the good agreement with the believed exact value 
ly = 3/4 and with the hyperscaling relation = 2A4 — 7, as soon as Nmin ^ 1000. 
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d = 2 SAW, Pure Power-Law Fit 



0.194 



0.192 



a 0.190 



0.188 



0.186 



10' 



10' 



10' 



Figure 2: Estimated exponent p from pure power-law fit to pivot-algorithm accep- 
tance fraction /, plotted versus Nmin- 
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d = 2 SAW Interpenetration Ratio 



0.68 



0.67 



^ 0.66 




0.65 



0.64 



10' 



10' 



10' 



10' 



N 



Figure 3: Interpenetration ratio \I/ versus iV, for 2-dimensional SAWs. Note that \I/ 
varies little for N ^ 100, and is constant within error bars for N ^ 1000. 
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0.592 



0.590 



d = 3 SAW, Pure Power-Law Fit 



0.588 



0.586 



0.584 



10^ 10^ 10^ 



Figure 4: Estimated exponent v from pure power- law fits to {R^y (□) and (i?^)^ 
(O), and estimated exponent (2A4 — 7)/3 from pure power-law fit to (T) (*), plotted 
versus Nmin- Note the very strong corrections to scaling, which lead to erroneous 
exponent estimates unless one takes Nmin ^ 10'*. 
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d = 3 SAW, Pure Power-Law Fit 



0.108 



0.106 



0.104 



10' 



10' 



10' 



Figure 5: Estimated exponent p from pure power-law fit to pivot-algorithm accep- 
tance fraction /, plotted versus Nmin- 
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0.592 



d = 3 SAW, Pure Power-Law Fit 



0.590 



^ 0.588 




0.586 



0.584 I — 
0.00 



0.05 

N 



0.10 



-0.5 



mm 



Figure 6: Estimated exponent v from pure power-law fits to {ReY ('-') ^i^d {RgY 
(O), and estimated exponent (2A4 — 7)/3 from pure power-law fit to (T) (*), plotted 
versus N~'^^ . Note the very roughly linear behavior, in agreement with the belief 
that Ai 0.5. 
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d=3 SAW, N,„,„=100, Level > 32% d=3 SAW, N„,„=150, Level > 33% 




0.585 0.5B6 0.5B7 0.5B8 0,589 0.585 0.586 0.587 0,588 0,589 



d=3 SAW, N„,„=300, Level > 32% d=3 SAW, N„,„=300, Level > 32% 




d=3 SAW, N„,„=400, Level > 32% d=3 SAW, N„,„=500, Level > 32% 




Figure 7: Pairs (u, A) for which the fit to the Ansatz ANp°'"^'' + BNp°'"^''-^ pro- 
duces a acceptable at the 32% significance level (i.e. one standard deviation). 
Observables are {RgY (^) ^-nd (T) (*). 
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d = 3 SAW Interpenetration Ratio 



0.28 



0.27 



^ 0.26 



0.25 



0.24 



10' 



10' 



10' 



10' 



N 



Figure 8: Interpenetration ratio \I/ versus iV, for 3-dimensional SAWs. Note that 
\I/ is a decreasing and convex function of iV, in flagrant disagreement with the 
prediction of the two-parameter renormahzation-group theory. 



83 



d = 3 SAW Interpenetration Ratio 



0.28 



0.27 



^ 0.26 



0.25 




0.24 



0.00 



0.02 



0.04 

^-0.56 



0.06 



0.08 



Figure 9: Interpenetration ratio \I/ versus N 
regression line is ^' = 0.24707 + 0.39312/iV°-^6 
plot. 



for 3-dimensional SAWs. The 
Note the excellent linearity of the 
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Figure 10: Wilson-de Gennes-type renormalization-group flow on the critical sur- 
face. H* (resp. Hq^) is the Gaussian (resp. good-solvent) flxed point. Ais (resp. 
Aiu) is the stable (resp. unstable) manifold of H*. Case la: Models in the good- 
solvent regime may have correction-to-scaling amplitudes that are either negative 
(P, -R) or positive (Q). Case II: The initial Hamiltonians Hn approach the stable 
manifold, while the low-energy effective Hamiltonians H!^^^ = TZ^Hn approach the 
unstable manifold. 
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d = 3 SAW, N = 70000. Bin size = 100000 



2x10' 



IxlO^ 



0x10' 



0x10^ 1x10^ 2x10^ 3x10^ 4x10^ 

Pivots 



d = 3 SAW, N = 70000. Bin size = 100000 



1.0x10" 



0.8x10" 



0.6x10" 



0.4x10" 



0.2x10" 




0.0x10" 

0x10*^ 



1x10*^ 2x10*^ 3x10® 4x10*^ 
Pivots 



Figure 11: Averages of (a) Rl and (b) T over bins of width 10^ iterations, for pivot 
algorithm on simple cubic lattice at iV = 70000 with parallel-rod start. 
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d=2 SAW. N = 1000, Autocorrelation Function 



d=2 SAW. N = 1000, Autocorrelation Function 




10 100 
t (pivots) 




10 100 
t (pivots) 



d=2 SAW. N = 1000. Autocorrelation Function 



^ 1 n-2 




S 10"^ ^ 



10 100 
t (pivots) 



;ure 12: Log-log plot of the sample autocorrelation function for (a) A = 

A = and (c) A = TKarp-Luby,R=20 for pivot algorithm on square lattice at 
= 1000. 
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